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INSPECTED 


-We  show^that  solutions  of  the  Cauchy  problem  for  systems  of  two 

conservation  laws  decay  in  the  supnorm  at  a  rate  that  depends  only  on  the 

1 

L/*  norm  of  the  initial  data.  This  implies  that  the  dissipation  due  to  the 

entropy  dominates  the  nonlinearities  in  the  problem  at  a  rate  depending  only 

V  A 

on  the  L r  norm  of  the  initial  data.^  Our  results  apply  to  any  BV  initial 
data  Uq  satisfying  Uq(±  ®)  =  0#  ahd  Sup{ug(»)}  <<  1.  The  problem  of 
decay  with  a  rate  independent  of  the  support  of  the  initial  data  is  central  to 
the  issue  of  continuous  dependence  in  systems  of  conservation  laws  because  of 
the  scale  invariance  of  the  equations .  Indeed,  our  result  implies  that  the 
constant  state  is^stable  with  respect  to  perturbations  in  L^^.  This  is  the 
first  st^bilirty  result  in  an  I<P  norm  for  systems  of  conservation  laws.  It 
is  <Jrucial  that  we  estimate  decay  in  the  supnorm  since  the  total  variation 

does  not  decay  at  a  rate  independent  of  the  support  of  the  initial  data. 

V 

^The  main  estimate  requires  an  analysis  of  approximate  characteristics  for 
its  proof.  A  general  framework  is  developed  for  the  study  of  approximate 
characteristics,  and  the  main  estimate  is  obtained  for  an  arbitrary  number  of 
equations.  '■ 
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DECAY  WITH  A  RATE  FOR  NONCOMPACTLY  SUPPORTED  SOLUTIONS  OF  CONSERVATION  LAWS 


Blake  Temple 


|1.  INTRODUCTION 


Consider  the  Cauchy  problem  for  a  system  of  n  conservation  laws 

ut  +  f(u>x  »  0  , 


(1) 


U(X,0)  =  UylX) 


where  u  =  (u1t...,un),  and  x  e  R,  t  e  R  •  We  study  decay  and  continuous  dependence  in 

solutions  of  ( 1 )  which  are  obtained  as  limits  of  approximate  solutions  generated  by  the 

random  choice  method  of  Glimm  [6 ] .  Thus  we  are  interested  in  solutions  that  take  values  in 

a  neighborhood  U  of  some  constant  state  u.  we  assume  that  df,  the  matrix  derivative 

of  f,  is  smooth,  has  real  and  distinct  eigenvalues  <  X^  <•••<  Xn  in  U,  and  that 

VX  *  R  >0  in  U  [9] .  Here  R„  denotes  the  unit  right  eigenvector  corresponding  to 
P  P  P 

eigenvalue  X^.  By  changing  the  frame  or  translating  the  flux  function  f  if  necessary, 
we  assume  without  loss  of  generality  that  u  =  0  and  X^  >  0,  p  = 

Let  u(x,t)  denote  a  weak  solution  of  (1)  which  is  a  limit  of  approximate  solutions 
generated  by  the  random  choice  method.  The  main  result  of  this  paper  is  the  following 
theorem  which  is  proved  in  the  case  n  «  2: 


Sjj§Qg§h  (1)  For  every  V  >  1  and  0  <  a  <  1  there  exists  constants  6  =  6(V)  <  1  and 
C(<J)  >  1  such  that,  if  uQ(*)  satisfies 

(2i  uQ(t  ">  =  0  , 

(3)  TV{Uq ( * )}  <  V  , 
and 

(4)  lu0(*)»s  <  6  , 

then 


1 

(5)  »uC,t)ls  <  ctoiilogi-^-A.,-  ]}  2+0 

L 

for  all  t  >  IUq(*)I  Here  constants  depend  only  on  f  and  their  arguments,  #  lg 
L 

denotes  supnorm  and  I  I  ^  denotes  L1  norm)  i.e. 

L 
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If  we  fix  the  initial  data  uQ ( • )  and  let  t  *  “,  then  (5)  gives  the  decay  of  the 
solution  u( • ,t)  in  the  supnorm  at  a  rate  independent  of  the  support  of  the  data.  Said 
differently,  (5)  verifies  that  the  dissipation  in  solutions  of  (1)  due  to  increasing 
entropy  overcomes  the  nonlinearities  in  the  problem  at  a  rate  depending  only  on  the  L^-norm 
of  the  initial  data  [cf  [9]].  If  we  fix  t  and  take  a  sequence  of  initial  data  tending  to 
aero  in  L1 ,  then  (5)  gives  a  rate  at  which  the  supnorm  at  time  t  tends  to  aero  with 
the  L1  norm  of  the  initial  data.  Because  the  values  of  u  at  time  t  have  a  bounded 
domain  of  dependence,  (5)  also  gives  a  rate  at  which  u(*,t)  tends  to  aero  in  as 

the  initial  data  tends  to  aero  in  This  is  the  first  continuous  dependence  result 

for  systems  in  the  norm  L1 . 

Other  decay  results  for  systems  have  bean  obtained  by  Glimm/Lax,  DiPerna  and  Liu  [4, 

5,  7,  10-14].  For  these  results  decay  is  obtained  by  means  of  estimates  for  the  decay  of 
the  total  variation.  In  the  case  of  nonperiodic  initial  data,  a  rate  of  decay  in  the  total 
variation  is  obtained  only  in  the  presence  of  compactly  supported  data,  and  the  rate 
depends  on  the  support  of  the  <bta.  It  is  crucial  in  (5)  that  we  estimate  the  decay  in  the 
supnorm  instead  of  the  total  variation  norm  because  simple  examples  show  that  the  total 
variation  does  not  decay  at  a  rate  that  depends  only  on  the  L1  norm  of  the  initial  data. 

Our  interest  in  the  L1  norm  in  (5)  stems  from  am  interest  in  the  problem  of 
stability,  by  which  we  meam  the  problem  of  the  continuous  dependence  of  solutions  on  the 
initial  data.  To  put  the  issue  of  stability  into  perspective,  we  make  the  following 
definitions:  we  say  that  solutions  of  (1)  are  strongly  stable  in  a  norm  I  I  if  there  is 
a  constant  C  >  0  such  that 

(6)  Ku(*,t)  -  v(»,t)l  <  Clu( • ,0)  -  v( • ,0) I  , 

for  all  weak  solutions  u  and  v.  We  say  that  solutions  of  (1)  are  weakly  stable  in  I  I 
with  a  rate  if 
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(7)  luC,t)  -  v( • #  t )  I  <  F(lu<*,0)  -  v(  •  #  0 )  I  • 

for  all  weak  solutions  u  and  v,  where  F  is  a  fixed  function  satisfying  lim  F(£)  »  0. 
It  is  well  known  that  solutions  to  scalar  conservation  laws  are  L  -contractive  [cf  8].  in 
the  language  above  this  says  that  solutions  of  a  scalar  conservation  law  are  strongly 
stable  in  L1  with  constant  C  =  1 .  Moreover,  by  studying  solutions  containing  a  single 
shock,  it  is  apparent  that  solutions  of  (1)  are  not  weakly  stable  in  the  supnorm  or 
in  the  total  variation  norm,  and  are  not  strongly  stable  in  any  lP,  p  >  1.  This  leaves 
L1  as  a  leading  candidate  for  studying  stability  in  systems  of  conservation  laws. 

Estimate  (5)  proves  that  the  constant  state  is  weakly  stable  in  Ljoc*  a  further 
comment,  in  [22]  it  is  proven  that  solutions  to  systems  are  not  L1 -contractive  relative  to 
a  constant  state  in  any  metric  that  is  compatible  with  the  u-space  topology.  This  directly 
implies  that  there  is  no  metric  D  and  constant  u  >  0  for  which  the  following  Gronwall 
type  estimate  holds  in  any  neighborhood  of  u: 

■jjT  D(u(x,t)  ,u)dx  <  oi  L>(u(x,0),u)dx  . 

Thus,  (5)  gives  the  stability  of  the  constant  state  in  in  a  regime  where  a  Gronwall 

inequality  fails  in  some  essential  way.  It  is  an  open  problem  whether  the  constant  state 
is  strongly  stable  in  LJoc* 

In  the  case  of  periodic  initial  data,  Theorem  (1)  holds  with  uQ( • )  replaced  by  the 
initial  data  in  one  period.  Thus  (5)  gives  a  decay  in  L°°  at  a  rate  independent  of  the 
period.  Again,  for  periodic  initial  data,  the  total  variation  does  not  decay  at  a  rate 
depending  only  on  the  L1  norm  of  the  data,  and  for  the  previous  decay  results  the  rate  of 
decay  in  the  total  variation  depends  on  the  length  of  the  period  in  the  initial  data.  Our 
methods  also  give  directly  that  periodic  data  decays  like  where  P  is  the  length  of 

a  period.  This  however,  is  not  sharp  in  light  of  the  t_1  decay  rate  obtained  by 
Glimm/Lax  [7]. 

We  now  indicate  the  proof  of  Theorem  (1).  Theorem  (1)  is  obtained  by  estimating  the 
decay  of  the  quadratic  functional  which  was  constructed  by  Glimm  in  [6].  Specif  real  Zy, 
let  h  denote  a  mesh  length  in  x,  and  let  uh(x,t)  denote  a  corresponding  approximate 
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solution  generated  from  initial  data  uQ(‘)  by  the  random  choice  method.  Roughly 
speaking,  the  values  of  u*1  at  time  t  are  obtained  by  approximating  the  actual  solution 
by  a  set  {y^}  of  simple  waves  each  of  which  moves  at  close  to  characteristic  speed.  The 
function  Q(t)  is  defined  by 

(8)  C(t)  =  I|Yt|  | Ym |  , 

where,  again  roughly,  the  sum  is  over  all  pairs  of  waves  at  time  t  that  will  interact  at 
some  later  time  due  to  differences  in  the  wave  speeds.  (In  the  words  of  Glimm,  summed  over 
all  "approaching”  waves.)  Here,  |y^|  in  (8)  denotes  the  strength  of  the  wave  y^  (for 
details  see  aections  (2)  and  (3)).  In  [6]  it  is  proved  that  Q  is  a  positive  decreasing 
function  of  time.  Heuristically ,  this  is  because  a  term  is  lost  from  the  sum  in  (8) 
whenever  two  waves  cross  each  other  in  the  xt-plane.  The  functional  Q  measures  the 
potential  for  interaction  of  waves,  but  contains  no  information  regarding  the  time  at  which 
interactions  will  occur.  Theorem  (t)  is  a  corollary  of  the  following  technical  lemma  which 
is  a  sharp  estimate  for  the  rate  at  which  Q  decreases  as  a  function  of  the  supnorm  and 
the  L1  norm.  For  this  lemma,  assume  that  n  ■  2  (or  that  there  exists  a  coordinate 
system  of  Riemann  invariants). 


LEMMA  (2A):  Let  Q  denote  the  quadratic  functional  associated  with  an  arbitrary 


approximate  solution  uh  which  is  generated  from  initial  data  ( • )  that  satisfies 


conditions  (2)  -  (4)  of  Theorem  (1).  Then  there  exists  a  constant  G  >  0  depending  only 
on  f  and  V  and  a  mesh  length  hg  «  hg(e,M)  such  that,  if 

(9) 


(10) 
and 

(11) 
then 

( 12) 


•»h(*'0)'s  >ii  ' 


I un ( • , 0 ) I  1  =  e  , 
L 


h  <  h. 


0  > 


g(0)  -  £> ( e ( gm ) 2 )  >  — 


(GM) 


.4. 


v ;  u1.1  r.1  w:  r 


In  words,  (12)  states  that  0  will  decrease  by  an  amount  on  the  order  of  the  supnorm 
squared  in  a  time  which  Is  on  the  order  of  the  L^-norm  divided  by  the  supnorm  squared.  For 
the  case  n  >  2  we  obtain  (12)  under  the  assumption  that  the  total  variation  of  uQ(*)  is 
small  [cf.  6] . 

LEMMA  (2B):  If  n  >  2,  then  there  exists  V  <<  1  such  that,  if  uQ(*)  satisfies  (2)  and 
(1)  of  Theorem  (1),  then  the  conclusions  of  Lemma  (2A)  hold. 

In  the  case  of  periodic  initial  data,  Lemmas  (2)  hold  with  uQ(*)  replaced  by  the 
initial  data  in  each  period. 

The  proof  of  Lemma  (2)  is  given  in  section  6.  (See  Theorem  6.3  for  a  detailed 
restatement  of  Lemma  2. )  The  proof  is  quite  technical  and  uses  the  theory  of  wave 
tracing.  The  theory  of  wave  tracing  was  developed  by  Liu  to  prove  that  the  random  choice 
method  converges  weakly  so  Long  as  the  sample  sequence  is  equidistributed  [15].  Wave 
tracing  is  a  method  of  keeping  track  of  left  and  right  states  on  approximate 
characteristics  [1,  4,  5,  ,  10,  15,  16].  Previous  decay  results  for  systems  use  the 

theory  of  approximate  characteristics,  but  rely  on  global  mechanisms  and  do  not  require 
keeping  track  of  left  and  right  states.  (It  is  important,  however,  to  recognize  that  in 
[5] ,  these  methods  are  localized,  and  decay  in  Q  is  used  to  control  decay  in  the  total 
variation  for  non-compactly  supported  data.  Of  course,  no  rate  can  be  obtained  for  decay 
in  the  total  variation.)  Here  we  develop  the  theory  of  wave  tracing  from  what  we  believe 
is  a  simpler  set  of  definitions  and  a  simpler  notation  than  has  been  previously  given.  The 
presentation  is  general,  and  essentially  self  contained.  Motivations  for  the  constructions 
can  be  found  in  [15,  16]. 

We  now  deduce  Theorem  (1)  from  Lemma  (2A)  using  the  basic  results  of  Glimm.  The 
remainder  of  this  paper  is  then  devoted  to  the  proof  of  Lemma  (2).  We  first  give  a  precise 
statement  of  the  results  in  [6].  (See  [23]  for  a  proof  of  the  supnorm  estimates.) 

LEMMA  (GL):  Assume  the  uQ  satisfies  the  conditions  (2),  (3)  and  (4)  of  Theorem  (1). 

Then  each  approximate  solution  u11  is  defined  for  every  h  >  0  and  t  >  0  arid,  moreover, 
there  exists  Gg  >  0  such  that 
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(13) 


TV{uh(*,t)J  <  aQ  TV{u0(*)}  , 

(14)  luhC,t)lg  <  GqIu0(.)Is  , 

(15)  luh(»,t2)  -  uh(» ,t1 )l  ,  <  GQ{h  +  |t2-t1 | >  , 

L 

and 

(16)  fl(t2  -  Q{ty)  <  0 

for  all  tj  <  tj*  (Proa  her*  on  out  we  uee  G0  to  denote  a  generic  constant  that  depends 
only  on  V  and  f • ) 

In  the  case  of  arbitrary  n,  the  results  in  [6)  are  that  (13),  (15)  and  (16)  hold  bo 

long  as  V  is  sufficiently  small-  Thu  reason  we  can  obtain  (5)  in  the  case  n  -  2  and 

not  n  >  2  is  that  we  use  (14),  a  result  that  is  not  known  for  n  >  2,  (For  n  >  2  one 

can  show  by  our  methods  that  thorn  exists  a  sequence  of  times  tj  ♦  »  for  which  (5) 
holds. ) 

So  assume  that  T.amms  (2A)  and  the  assumptions  of  Theorem  (1)  hold.  Let  M  >  1  be 
given,  we  estimate  the  time  at  which  lu^(»,t)lg  <  for  h  <  h0-  Let  Gq  >  1  be  large 
enough  so  that  Q(0)  <  Gq.  Set 

(17)  G)  3  GqG 
and  let 

(18)  N  *  [(G1M)2]  +  1 

where  [  ]  denotes  "greatest  integer  in".  Define  the  times  t ^ , . . . , tN  between  which  Q 
decreases  by  an  amount  (G ■,»*) ~2  as  follows) 

(19)  t0  *  0  , 

tn+1  i  «up{t  >  tn  ,  C(tnf)  -  C(t)  <  — '—s)  • 

Define 

(20)  e  5  luh(«,t  +  )l  . 

n  n  i 

L 

Tint  m  <  N  lie  that  smallest  integer  for  which 

<*1>  Vtl  "  >  Em(GiM)2 

if  such  an  integer  exists-  Otherwise  let  m  «  N.  Now  if  t^  <  •»,  then  Lemma  (2A)  implies 
that 

<22>  ,uh<*'t»+,,S  ‘  G^M  • 
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To  see  this,  note  that  the  contrapoaitive  of  Lemma  (2A)  states  that  if 

Q(0)  -  fc(e(GM)2)  <  — —•  ,  then  luh(‘,0)l_  <  -jj.  if  t_  <  then  this  applies  with 

h  (GMih  ”  * 

u(«,0)  replaced  by  u  (*«ta+)  and  M  replaced  by  G0M.  (Note  that  if  m  »  N  and 

tN  <  •»,  then  (21)  holds  in  this  case  because  Q  can  incur  no  more  than  h  decreases  of 

magnitude  (GM)“2. )  Thus  if  <  «,  then  by  (14) 

■uV.tMg  4  i 

for  all  t  >  tg,  as  stated. 

It  remains  only  to  estimate  t^.  We  show  that 

(G  M)2 

(23)  tffl  <  (G2M)  e0 

where  G2  *  2G0G.|.  Without  loss  of  generality,  assume  that  t1  >  e  and  t1  >  h  so  that 
(15)  gives 

124)  en  <  G0tn  , 

where  again  we  take  G0  to  be  generic.  Thus  if  n  <  m, 

(25)  t„+1  -  t„  <  <  G0tn { G jM ) 2  , 

so  that 

(26)  t„+,  <  {1  +  G^(GjM)2}tn  . 

By  induction  this  implies  that 

tm  «  (1  +  G0<G1M)2}'B"1t1 

<  (1  +  G0(G1M)2)m"1E0(G1M)2 

<  (1  +  G0(G1M)2}“e0 

(27)  <  {1  +  G0(G1M)2}NeQ 

2  (G1M)2 

<  {1  +  Gq(G1M)  }  Eg 

<G  M)2 

<  (G2M)  Cg 

where  G2  =  2G0Gf  Thus  we  have  that 

(28)  luh(*,t)«s  < 
so  long  as 
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In  particular,  let  t  *  (G2M) 
(29) 


(G2M)" 


(G  M) 

t  >  (G2M)  e0 


and  choose  0  <  a  <  1.  Then 

2+a 


log[— J  -  (G  M)2log(G  M)  <  (C(O)M) 
60 


for  some  C(o)  >  0.  Thus 


C<C){log(i-j}  2+0  >  1  >  luh(-,t)I  . 
e0  ” 

Now  if  u  is  a  limit  of  approximate  solutions  uh  as  h  +  0,  then  it  follows  that 

luh(*,0)l  ♦  lu  (•>»  and  we  conclude  that 

L  L 


(31)  luC,t)ls  <  C(o)(log[  lu^(^) ,  Jl 

This  completes  the  proof  of  Theorem  ( 1 ) . 


1 

2+0 


-°-  '’J 


In  the  case  of  periodic  data,  the  estimate  (31)  is  obtained  in  the  same  manner  by 
replacing  u0  by  the  initial  data  in  one  period.  In  this  case  C(o)  is  independent  of 
the  period.  For  periodic  data  we  can  also  use  Lemma  (2)  to  obtain  a  rate  of  decay  which 
depends  on  the  period,  since  for  periodic  data. 


(32) 


lu  ( • ,t) I  1  <  G0P 
L 


where  P  is  the  length  of  one  period.  In  this  case  we  can  use  (32)  instead  of  (24)  and 
argue  as  follows:  if  uh(*,t)  >  g-jj,  then  Lemma  (2)  gives 

®  .  i 

Q(t)  -  Q(t  +  G  P(G  M) i  )  >  - —  . 

(G^)2 

2 

Since  Q(0)  <  G0,  we  see  that  Q  can  incur  only  G0(GiM)‘  decreases  of  magnitude 


<GtM)' 


before  Q(t)  < 


r»  in  which  case 


(G1M) 


lubc.tXg  <  ^  • 

Thus  lu^t'.tllg  must  be  smaller  than  —  before  time  T  =  GqP(G1M)2Gq(G1M)2  =  CPM4. 
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This  gives  that 


for  all  t  >  CPM  ,  in  which  case 


"uh(*'t)as  < 


Huh(.,t)lc  <  C(|)“1/4  . 


We  note  that  (33)  is  not  sharp  (t-1  is  sharp,  cf.  Glimm/Lax  [7]).  This  might  be  expected 

since  we  are  not  invoking  global  mechanisms  of  decay  as  in  [7]. 

The  remainder  of  this  paper  is  devoted  to  the  proof  of  lemmas  (2)  and  (3).  Before 

embarking  on  the  proof,  we  briefly  discuss  the  idea  behind  it.  The  idea  is  that,  since 

Uq(±  °°)  =  0,  if  *uo^*^*s  =  m"  an£*  "V”  i  =  £'  then  there  must  be  a  "spike"  in  the 

L 

initial  data  of  height  on  the  order  of  —  in  | u |  and  width  in  x  on  the  order  of  £M. 

M 

For  example,  if  |uB(x)|  >  for  all  x  e  [xA,xB]  ,  then  e  =  II  uQ  ( •  )  M  1  >  |  x0  -  xA| 

L 

which  implies 

|xB  -  xA|  <  2Me  . 


Thus  consider  the  case 


XA  <  x  <  *B 


otherwise 


where  | vi |  *  — ,  xB  -  xA  =  eM.  This  data  resolves  into  four  simple  waves  associated  with 


the  Riemann  problems  [0,u]  and  [u,0]  [cf.  §2,  [9)).  Label  these  waves  o^,  o^,  S1  and 


as  in  figure  1. 


Figure  1 
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has  been  defined  in  such  a  way  that  properties  (4.1)  -  (4.3)  are  satisfied. 

Mp(J+1)»  r^tJ+l)  in  terms  of  Mp(J)/  rp(J)  and  see  that  properties  (4.1)  - 

for  Mp(J+1)»  Tp(j+1),  respectively. 

To  define  M„(J+1)  and  T  (J+1),  we  define  the  sets 
P  p 

M  £(  J+1 )  =  U  eM  (J+1)  :  t(J)  =  i}  , 

P  p 

and 

{'(jj}  A  » 

JleMp(J+i ) 

for  each  i  e  a.  The  sets  MfJ+1)  and  T  (J+1)  are  then  defined  by 

P  p  * 

M  (J+1)  =  (M  (J>  n  M(J-D)  u  (u  Mi(J+D)  , 

P  P  ^  P 

rp(j+i)  =  (Yg  :  t  e  Mp( J+i ) }  . 

So  fix  i  e  Z.  Define 

(4.12)  L.  =  (ted  (J)  :  i(J)  =  i,  sgnU)  ft  Bgn(yP,)}  . 

P  *J 

If  «  0,  then  define 

il£(J+1)  =  L*  > 

P  * 

r£(j+u  =  (ij  :  t  e  L*>  • 

So  assume  ft  0.  let  {iQ  ' " '  *  '*a-1  ^  —  ^p(J)  denote  the  indices  in 

satisfying 

*k<J>  -  i  » 


and  ordered  so  that 


sgnU^) 


Vi  4  V  1  4  k  <  a_1' 


»L  E  L(YiJ>  ' 


uRHR(Yf.)  , 


k-1 


sgn(Y?j)  , 

in  the  sense  of  Property  (4.1) 


oQ  =  °,  ok  =  l  yt  (tj-),  k  =  0, . . . ,a  , 
s=0  s 

uk  E  T<V  V  • 


We  define 
(4.3)  hold 


Wp(J) 


.  Define 
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Property  (4.3)  expresses  the  fact  that  characteristics  trace  nonzero  elementary  waves 
of  a  given  family  and  sign: 

PROPERTY  (4.3):  If  l  e  Mp(J),  then  the  signed  strength  Yj(t)  of  the  characteristic 

0  1 

Y^  is  constant  and  nonzero  in  (t^t^).  We  write 


(4.11) 


For  convenience,  we  set 


sgn{y^(t) }  =  sgn( i)  . 


YA(t)  -  0 


for  t  4  [t°,tj). 

We  now  define  Mp(j)  and  r^tJ)  by  induction  on  J.  We  simultaneously  verify 
properties  (4.1)  -  (4.3)  which  are  assumed  in  the  induction  step. 

First  assume  J  =  1.  Define 


Mp(1)  “  Uj  :  ^  t  0,  p  -  1,2,  q  =  L,R} 


where  is  defined  by 


»^(0)  -  i  , 


*?(1)  -  ( 

^  i+1 


if  q  =  L 
if  q  -  R 


Then  for  0  (  t  <  k  i  tt,  and  /  0,  define 


U^q(t)  "  * 

i 


-*(t,  -R(Y«)  • 


It  is  easy  to  verify  that  properties  1,  2  and  3  hold  for  Mp(D  and  Tp(1). 

Now  assume  for  induction  that  Mp(J)  £_M(J)  has  been  defined,  and  that  for  every 


l  e  Mp ( j ) , 


Yjtt]  =  (u^(t)  ,u*(t>  ) 
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CR  =  (so  that,  e.g.  uR  =  ^p(°R!UL^*  Then  for  tj  <  t  <  tj+i» 


u*  (t)  “  uh  ' 
0 


(t>  *  UR  ' 
a 


UJ  ,(t)  =  U^(t)  =  T<  W  ' 

k- 1  k 

where  e  (0,oR)  end  |ok+1|  >  |ok|,  1  <  k  <  a-1.  We  define 

(4.6)  Ta  (t)  3  ak+1  -  ok  . 

k 

(i.e.,  round  brackets  around  t  to  distinguish  it  from  Y^ft) )  to  be  the  signed  strength 
of  the  characteristic  y4  at  time  t  e  [tj,tj+1),  1  <  k  <  a-1.  Moreover,  for  a  e  (0,1) 
and  t  e  [tj,tj+1),  we  define 

(4.7)  u«  (t)  =  T(ak  +  <J«Jk+1  ~  \)>  uL) 

k 

for  1  <  k  <  a-1.  Note  that  if  Property  (4.2)  holds,  then  (4.6)  and  (4.7)  define  y^tt) 

and  u°(t)  for  all  l  e  Mn(J),  t  <  t,.  We  use  (4.7)  to  define  two  characteristics  y  T 
*  p  l 

o 

and  y  D  corresponding  to  each  a  e  (0,1)  and  each  i  e  M„(J)  satisfying  t(J)  *  i,  as 
follows : 


t(j)  if  j  <  J-1 


*’(j)  -  <  i  if  j=J,  q  =  L 

V.  i+1  if  j  =  J,  q  *  R 

Y  L [t]  =  (u£(t),u°(t)  )  , 


(4.10) 


Y  R[t]  =  <u°(t),u*(t)>  . 

*0 


(I.e.,  o  determines  a  splitting  of  the  characteristic  y.  into  a  characteristic  y  of 

* 

o 


strength  a  Yj(t)  and  a  characteristic  y  R  of  strength  ( 1 —a >y ^ ( t ) ) . 
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M(J)  -  u  M  (J)  , 

P-1 

the  RHS  being  a  disjoint  union.  Corresponding  to  each  i  e  Mp(J)  is  the  characteristic 

■  <ut'uf>  e  rp(J)  ' 

each  entry  being  a  function  of  t  for  t®  <  t  <  <  tj.  we  let  y^tt]  =  (u£[t] ,u* [t] ) 

L  R 

denote  the  value  of  y^  at  time  t.  The  functions  u^[tj,  u^[t]  (the  left  and  right 
states  of  the  characterisitcs  at  time  t)  are  constant  on  intervals  tj  <  t  <  tj+1 

for  j®  <  j  <  jj.  r ( J )  is  the  disjoint  union 


r<J)  -  U  T  (J)  . 

P-1  P 

For  convenience  we  set  Y^It]  «  0  for  t  ft  [t° ,  t  J ) . 

Before  defining  M( J )  and  T(J)  precisely,  we  first  list  three  properties 
(properties  (4.1),  (4.2)  and  (4.3))  which  the  characteristics  satisfy.  Then  for 
p  -  1,...,n  we  simultaneously  define  Mp(J),  T^(J),  and  verify  properties  (4.1)  -  (4.3) 
by  induction  on  J. 

Property  (4.1)  states  that  each  set  Mp(J),  p  »  1,...,n  is  partially  ordered,  and 
expresses  the  fact  that  characteristics  of  the  same  family  never  cross t 


PROPERTY  (4.1):  If  e  Mp(J),  then  for  every 

(4.4)  j  €  IjJ  ,j!  1  n  lj®  1  , 

*1  »2  l2 

we  have  that  either  I,(j)  <  i,(j)  or  t,(j)  <  t.(j).  If  (j?  ,j]  1  n  lj°  ,j]  ]  ¥  9  and 
1  i  1  l1  l2  *2  1 
t^(j)  <  i2(j)  for  J  satisfying  (4.4),  then  we  say  that  f,  <  l2  (or  lies  to  the 

left  of  fj). 

Property  (4.2)  states  that  the  left  and  right  states  of  the  p-characteristics 
l  e  l^,(J)  that  satisfy  I(j)  “  i,  j  <  J,  partition  the  p-wave  curve  Tp(*,L(y^)) 
between  L(y*j\)  and  R<^ij  > : 

PROPERTY  (4.2):  Let  i  e  Z,  j  8  [O.J-1]  be  such  that  y*\  f  0.  Let  UQ , . .  . ,  t^}  be 
the  set  of  p-characteristics  t  e  MpU)  satisfying  i(j)  »  i,  and  ordered  so  that 
*k  *  ^k+1  the  sense  of  Property  (4.1).  Let  uL  =  L11^j)»  uR  =  R(Y^j)  and  ^et 
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•  •J*. 


characteristics  in  detail. 


where  j°  <  j1  are  in  the  set  {0,1,..., j}.  Let 

(4.1)  S(  j°,  j 1 )  *  U  :  M(  j°,  j 1 )  +  Z  :  4(  j)  -  4(j-1)  e  {0,1}}  . 

For  each  element  of  1  e  M(j®,j^),  define  the  function 

(4.2)  Xj  :  [t  Q,t  11  ♦  R 

3  j 

as  follows: 

X4(V  *  x4(j)  '  j°  «  5  <  31  . 

xt(t)  *  +  >  “  xg(tj  )]h,  t^  <  t  <  tj+1  . 

Thus  each  element  4  e  M(j  >j  )  corresponds  to  a  continuous,  timelike,  piecewise  linear 
curve  in  the  xt-plane  given  by  the  graph  of  x^.  Note  that  the  graph  of  x^  is  defined  in 
[t  0,t  .,]  and  connects  successive  mesh  points*  i.e.,  the  slope  of  the  curve  is  either  0 
or  If  j,  t  is  outside  the  domain  of  (,  i(,  then  we  write  Jt{  j )  «  0,  x^tt)  “I i, 

respectively . 


(4.3)  Sf(J)  -  U  mj0,}1)  • 

0<j°<j 

Then  for  4  e  M(J)#  define  jj»  to  be  the  positive  integers  such  that 


4  e  tf(j°j!)  . 


0  1 

Let  =  t  Q  ,  tt  =  t  ^ 


We  presently  define  the  set  M(J)  of  indices  for  the  characteristics,  as  well  as  the 
set  T(J)  of  characteristics,  by  induction  on  J.  For  p  «  1,...,n,  the  set  of  indices 


for  the  p-characteristics  is  a  set 


Mp( J)  c  Sf(J)  , 


§4.  DEFINITION  OF  APPROXIMATE  CHARACTERISTICS 

In  this  section  we  define  the  set  of  approximate  characteristics  (heretofor  referred 
to  as  "characteristics” )  associated  with  a  given  approximate  solution  u*1  and  time  level 
Tj.  In  the  next  section,  we  study  properties  of  the  characteristics.  The  procedure  is  as 
follows.  He  first  define  an  index  set  M(J)  for  the  timelike  piecewise  linear  curves  that 
connect  successive  mesh  points  (xiftj)  in  the  xt-plane  [cf.  16].  A  subset  M(j) 
corresponds  to  the  set  of  characteristics.  We  call  this  the  index  set  for  the  character¬ 
istics.  Each  element  t  @  M(J)  gives  the  position  of  an  "elementary  wave"  y^  at 
different  time  levels.  The  piecewise  linear  curves  that  are  undefined  at  t  «  tj 
correspond  to  elementary  waves  cancelled,  and  those  undefined  at  t  *  0  correspond  to 
elementary  waves  created  by  nonlinearities.  He  let  H  M(J)  denote  the  index  set  for  all 
such  characteristics  [cf .  15 ] • 

He  define  the  elementary  wave  Y4  associated  with  t  e  )|(J)  by  assigning  a  left 
L  R 

state  [t]  and  right  state  u^ [t]  to  each  time  level  that  intersects  the  piecewise 

linear  curve  defined  by  i.  We  define  y.[t]  =  (u^[t]  ,u*[t]  ) ,  we  call  T(J)  -  U  {y.} 

11  JteM(J)  1 

the  set  characteristics  defined  for  time  tj.  The  assignment  of  states  to  characteristic 

curves  is  done  as  follows:  We  first  state  three  properties  that  the  assignment  should 

satisfy,  and  then  we  assume  the  properties  to  hold  in  order  to  define  the  characteristics 

at  the  induction  step.  Thus  the  characteristics  are  defined  and  the  properties  are 

verified  simultaneously  at  the  induction  step. 

The  characteristics  determine  a  partitioning  of  the  waves  in  uh  appearing  before 
time  tj.  It  is  important  to  estimate  the  "fineness"  of  the  partition.  For  Liu  [15],  the 
fineness  is  built  into  the  procedure  by  an  initial  partitioning  of  the  waves.  The  cost  of 
having  M(J)  as  an  index  set  is  that  we  must  estimate  the  fineness  of  the  partition  as  a 
function  of  J.  This  together  with  an  estimate  for  the  speed  of  characteristic  curves  is 
given  in  Theorem  (5.12).  The  remainder  of  section  5  is  essentially  devoted  to  obtaining 
estimates  for  the  change  |u^[t]  -  u^[0]|,  q  -  L,R,  in  terms  of  changes  in  Q  and  in 
terms  of  strengths  of  elementary  waves  assigned  to  characteristic  curves  that  cross  the 
characteristic  curve  of  l  in  time  [0,tj] .  We  now  proceed  with  the  definition  of  the 
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then  uh  is  well  defined,  takes  values  in  U  for  all  time,  and  moreover 
(3.18)  Vj  <  G0V„  , 

(3-19)  Gj1  1  Di;)  <  fi(tr)  -  Q(tj+)  <  Gq  l  Dij 

and 

(3.20)  luh(«,t)  -  uhC,s)l  1  <  GQ{h  +  |t-s|} 

L 

for  all  j,  s,  t  >  0. 

LEMMA  (GL3).  Assume  that  Uq(±«“)  »  0  and  that  there  exists  a  coordinate  system  of  Riemann 
invariants  (eg  n  *  2).  Often  for  every  V  >  0  there  exists  M  >  0  such  that  if 

(3.21)  V0  <  V  , 
and 

(3.22)  lu£«s  <  ^  . 

then  uh  is  well  defined  for  all  t  >  0,  takes  values  in  U,  and  (3.18)  -  (3.20)  hold 
together  with 

(3.23)  luh(*,t)1s  <  Ggluh( • ,0 ) lg  . 

(See  (23)  for  a  proof  of  Lemma  (GL3 ) . ) 
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V  . 


Two  waves  Y^j  and  y®ij  are  said  to  approach  at  time  tj  if  one  of  two  conditions 
holds  [6]  : 

p  <  p'  and  i  >  i*  , 


(3.10A) 

or 

(3.10B) 


p  «  p'  and  either  Y^j<  0  or  Yi'j  <  ° 


For  tj  <  t  <  tj+1,  define  [4] 

(3.11)  Q(t)  5  JlY^jl  lY^'.jl 

where  the  sum  is  over  all  pairs  of  waves  Y^j  and  Y^,j  that  approach  at  time  tj. 

Define 

(3.12)  °ij  =  Z|Y*?1>j_1l  lY^j.J 

where  the  sum  is  over  all  pairs  of  approaching  waves  that  enter  the  diamond  i  and 
finally  define  the  cancellation 

<=•'31  <?,  ,  i  tlYf!,,,.,!  .  U*  „,l  -  It*,,.,  *  Y^.,1  . 

Then  c£j  measures  the  amount  of  p-wave  cancelled  from  both  Y^*1jj_1  and  Yi^j-1  at  th® 
interaction  in  A 


ij 

The  following  lemmas  are  due  to  Glimm  [4} .  Let 

(3.14)  Vj  S  l  Wf  j  I 

i»P 

estimate  the  total  variation  of  uh(»,t)  for  t^  <  t  <  tj+i* 


LEMMA  ( GL 1 ) .  If  un  U0,  then  for  all  i,  j. 


I Yij  "Ti-1,j-1  ~  Yir  j-1 1  *  G0Dij 


(3.14) 

This  immediately  implies 

(3.16)  |Y?j|  -  lY^I  +  lY^j.,1  -  2^  +  0(1)0^ 

with  1 0(  1 )  |  <  Gq.  (Again,  Gg  denotes  a  generic  constant.) 


LEMMA  (GL2 ) .  There  exists  a  constant  V  >  0  and  a  number  Gg  >  0  such  that  if 
(3.17)  V0  <  V  , 
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(3.7) 


L(Y li  »  *  L<YiJ» 


R(YfJ  )  -  uh(x1  ♦  aj+1h.  tj+,->  -  L(Y£) 


R(y£R  )  -  R<Y^) 


More  generally.  If  (3.6)  does  not  hold,  then  define 


yPj  if  >  s(YPj) 


0  otherwise 


•  Yf*  <  0 


r?,  if  a,  <  XJL(Y?.))  ,  Y?,  >  0 


^  if  <  .(Y^)  ,Y^<0 


0  otherwise 


[Here,  for  example,  Y^  is  defined  to  be  zero  if  the  wave  lies  to  the  right  of  the  sample 
point  +  ajh,  and  Y^j  if  it  lies  to  the  left] . 

By  construction,  the  waves  Y^j  solve  the  Riemann  problem  for  »  uh(xi_1  + 

“j-lh,  tj)  =  L(Y^j>,  ur  ”  uh*xi  +  aj-th'  tj)  =  that  iB  Pos*d  at  (xA,tj)  in  the 

approximate  solution  uNx.t).  Because  we  assume  that  all  wave  speeds  are  positive,  the 

n  pR  pL 

waves  Y*  ^  are  formed  by  the  interaction  of  the  waves  Y^_1  j_1  with  the  waves  Y  ^ 

at  time  tj_j>  P  =  1,2.  To  emphasize  this,  we  let  6^  denote  the  diamond  of  interactions 

centered  at  (x^tj)  [6,  16]  *  i.e. ,  consisting  of  vertex  points  (x.^  ♦  aj.^h,  tj), 

1  1  P 

(xt  +  aj_.,h,  tj),  (*i#tj  -  y  k)  and  (x^tj  +  fc).  We  say  that  the  waves  Y^j  which 

cross  the  upper  wedge  of  the  diamond  are  formed  due  to  the  interaction  of  the  waves 

Y^i  j.j  and  Y^Lj_1  which  cross  the  lower  wedge  of  the  diamond.  We  call  YRj  a  wave 

pR  pL 

that  leaves  the  diamond  A. .,  and  we  call  the  nonzero  waves  among  yT  .  .  ,  and  y,  .  . 

13  i“i»3_1  i»3”' 

the  waves  that  enter  the  diamond  . 
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}3.  THE  RANDOM  CHOICE  METHOD  APPROXIMATES 


! 

| 

I 


In  this  section  we  define  the  approximate  solutions  of  ( 1 )  generated  by  the  random 
choice  method  of  Glims  [6],  and  develop  notation  required  for  the  subsequent  sections . 

Let  h  be  a  mesh  length  in  x(  and  let 

(3.1)  k  -  Ch 

be  the  corresponding  mesh  length  in  t,  C  >  Sup  {|X  ( u >  J } .  For  i,  j  e  2,  j  >  0,  let 

ueu,  n 

=  ih,  tj  s  jk.  Let  a  be  a  sample  sequence,  a  s  0  <  »j  <  For  given 

initial  data  uQ ( • )  c  u,  define  the  random  choice  method  approximate  solution 
uh(x,t)  =  uh(x,t »£)  by  induction  on  j  as  follows:  First,  for  xA  <  x  <  xi+1,  define 

(3.2)  uh(x,0)  i  uj(x)  -  uQ(xi  +  -|)  . 

Next,  assume  for  induction  that  u^(x,t)  has  been  defined  for  t  <  tj.  Define 

(3.3)  uh(x,tj)  =  uNxj+ajhjtj-)  , 

and  for  tj  <  t  <  tj+1,  define  uh( x,t)  to  be  the  solution  of  the  Riemann  problem  posed 
in  (3.3)  at  time  tj.  By  (3.1),  uh  is  well  defined  so  long  bb  uh(x,tj)  £  U  for  all 

V 

Let  uh  be  any  approximate  solution  that  is  well  defined  by  the  above  procedure.  Let 
Y^j  denote  the  name  as  well  as  the  signed  strength  of  the  p-wave  that  appears  in  the 
solution  of  the  Riemann  problem  that  is  posed  at  (x^tj)  in  the  approximate  solution 
uh.  Define 

L«1J»  2 

R(Yij)  H 

If  Y^j  is  a  shock  wave,  define 

(3.5)  s(Y^j)  2  speed  of  wave  Y^j 

If  Y^j  is  a  rarefaction  wave,  then  the  wave  is  "split"  at  time  fcj+i  if 

(3.6)  X  (L(Y?  , ) )  |  <  a.  <  X  (R(Y?j  ) )  £  « 

P  ij  h  j  p  ij  h 

In  this  case  define  Y?j  and  Y^j  by 


left  state  of  the  wave  Y^j  , 
right  state  of  the  wave  Y^j  • 
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•  -  %  • . 
.-i'-AV. 


by  the  condition 


uk+1  e  MV'  1  «  k  4  n  • 

We  call  Tp  [cf  [11)  a  wave  strength  measure  for  the  p-simple  waves  if,  for  each 

UL  e  °1«  Tp(o/uL)  is  *  parameterization  of  in  Uj,  a  e  R,  and  moreover 

Tp  e  C  with  bounded  third  derivatives  with  respect  to  both  arguments,  Tp(0;uL)  =*  uL  and 

3X  (T  ( c ,u  ) ) 

(2.3)  -E-f; - >  o  • 

If  LIT**)  ”  u^  and  R(YP)  “  Tptoju^),  then  we  define  the  signed  strength  of  yP  to  be 
o.  We  let  Yp  =  o  so  that  "yP"  denotes  both  the  name  as  well  as  the  signed  strength  of 
the  wave.  By  (2.3),  yp  <  0  for  shocks  and  Yp  >  0  for  rarefaction  waves. 

In  the  following  three  sections  we  study  the  random  choice  method  approximates  and 
associated  approximate  characteristics  using  arbitrary  T  to  define  wave  strength.  In 
section  b  we  define  T  by  means  of  a  best  approximation  to  a  coordinate  system  of'Riemann 
invariants.  The  following  lemma  is  a  direct  consequence  of  the  fact  that  the  speed  s  of 
a  shock  is  equal  to  the  average  characteristic  speed  to  within  termB  that  are  quadratic  in 
the  strength  of  the  shock  [9,  11,  17): 

LEMMA  (2.1).  Assume  that  uR  «  Tp(<JiuL>,  o  <  0.  Then  the  speed  s  of  the  corresponding 
shock  wave  satisfies 

u  +  u_  , 

s  -  Xp( - 5 - )  +  0(1)0 
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J2.  THE  RIEMANN  PROBLEM 

The  Riemann  problem  Is  the  initial  value  problem  ( 1 )  where  the  initial  data  has  the 

form 


/  ^ 

for 

X 

A 

O 

u0(x)  -  < 

K 

for 

x  >  0 

He  assume  that  the  Riemann  problem  is  uniquely  solvable  by  the  method  of  Lax  [9]  for  all 
uL  and  uR  in  a  neighborhood  U  of  u  ■  0.  In  particular,  assume  that  all  states  that 
appear  in  the  solution  of  a  Riemann  problem  posed  in  U  lie  in  a  set  2.  Assume  that 
Ap(u)  <  Ip, (v)  for  all  u,  v  e  1  <  p  <  p'  <  n,  and  that  (1)  is  genuinely  nonlinear  in 
all  characteristic  fields  throughout  U  1 . 

Let  Rp(u^)  denote  the  integral  curve  of  Rp  that  contains  the  point  uL.  Let 
Rp(uL)  denote  the  p-rarefaction  curve  associated  with  the  point  uL»  i.e.,  that  portion 
of  Rp(uL)  for  which  Xp(u)  >  X^(uL) .  By  [91  there  exists  a  unique  curve  Sp(uL)  that 
makes  C2  piecewise  contact  with  Rp  ( uL )  at  the  point  uL,  and  such  that  for  each 

u  e  Sp(uL)  there  is  a  scalar  s  such  that 
(2.1)  s [u]  -  [f]  , 

where  [u]  S  u  -  uL,  [f]  *  f(u)  -  f(uL).  Statement  (2.1)  is  the  Rankine  Hugoniot  jump 

condition,  and  we  say  that  Sp(uL)  is  in  the  Hugoniot  locus  of  uL.  Let  S  ~ ( uL )  denote 

the  p-shock  curve  associated  with  uL»  i.e.,  that  portion  of  Sp(uL)  for  which 

Xp(u)  <  Xp(uL) .  He  assume  that  Xp  is  monotone  on  both  Sp(uL)  and  Rp( uL) ,  so  that  the 

curve 

(2-2)  Vp(uL)  =  Sp(uL)  u  Rp(uL) 

is  a  C  piecewise  CJ  curve  throughout  U^.  For  uR  e  yp (uL) ,  the  Riemann  problem  is 
solved  by  a  p-simple  wave:  a  shock  wave  if  Xp(uR)  <  X^ (u^ )  and  a  rarefaction  wave  if 
X  (u„)  >  X  (u  )  [9,  17].  He  let  YP  denote  any  p-simple  wave.  For  given  yp,  let 

P  K  p  u 

I*(YP)  denote  the  left  state  of  YP  and  let  R(YP)  denote  the  right  state.  For  arbitrary 
states  uL,  uR  e  U ,  the  Riemann  problem  is  solved  uniquely  by  n  simple  waves  Y  \  . . •  ,Yn 
which  are  separated  in  the  xt-plane  (going  from  left  to  right)  by  the  constant  states 

P 

U1  2  uL'’**,un+1  E  UR*  T^le  states  u^,  and  hence  the  waves  Y  »  are  uniquely  determined 
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It  takes  the  longest  time  for  waves  to  Interact  if  a2  *  0  or  6  ^  —  0.  Assume  that 
<*2  «  0  and  to  be  specific,  assume  that  is  a  shock  wave.  In  this  case,  must  be  a 

rarefaction  wave  because  u  -  0  is  both  the  left  most  and  right  uiost  state  in  the  problem. 
(Here  we  use  the  assumption  that  the  waves  are  weak.)  Now  we  can  estimate  the  time  at 

which  fi  must  decrease  by  order  i.e.,  the  time  it  takes  to  interact  with  8 1 . 

M  1 
But  the  differences  in  speeds  between  and  8 1  at  time  zero  is  on  the  order  of  — ,  so 

interaction  must  occur  within  a  time  on  the  order  of  T  where  T  satisfies 

l*B  -  xa(  -  ^  *  “  0  * 


T  -  m|xb  -  XA|  -  M2E  . 

Since  two  waves  a1  and  8 1  of  strength  on  the  order  of  ^  interact  before  time  T,  by 

(8)  we  expect  Q  to  decrease  by  an  amount  on  the  order  -1—  in  time  T.  This  is 

1  M 

fi<0)  -  Q(T)  -  0(1)  , 

M 

which  is  (12). 

We  implement  the  above  idea  as  follows  >  given  data  u. ( * ) ,  we  locate  p-waves  a  and 


8  whose  speeds  at  time  t  ■  0  differ  by  0(1)  —  and  whose  x  distance  apart  at 

M 

t  «  0  is  0(1)  eM  as  above.  We  identify  these  waves  at  a  later  time  by  means  of 

approximate  characteristics.  We  then  assume  for  contradiction  that  Q  does  not  decrease 
1  2 

by  0(1)  —  in  time  0(1)  EM  .  A  consequence  of  this  is  that  a  and  8  may  be  chosen  so 

M  2 
that  the  corresponding  characteristics  do  not  intersect  before  time  0(1)  EM  .  We  finally 


derive  a  contradiction  by  estimating  that  since  the  decrease  in  Q  is  small,  the  speed  of 

the  characteristics  a  and  8  agree  with  the  speeds  at  time  zero  sufficiently  to 

2 

guarantee  that  they  intersect  before  time  0( 1 )  EM  .  By  this  contradiction  we  can  conclude 
the  proof  of  lemmas  (2)  and  (3).  In  the  above  analysis  we  must  keep  track  of  left  and 
right  states  on  approximate  characteristics.  This  is  essentially  the  wave  tracing  idea  of 
Uu  [15J. 

In  sections  2  and  3  we  review  the  Riemann  problem  and  the  random  choice  method  and 
establish  notation.  In  section  4  we  define  approximate  characteristics.  In  section  5  we 
establish  properties  of  approximate  characteristics.  (This  is  done  in  a  general  setting.) 


‘vlv-l v? I- ••5 


-.-y 


V-»V  > 


4n(5)  S 


0 

if 

j  <  J  , 

i 

if 

j  -  J  > 

i 

if 

j  =  J+1. 

o  <  a, 
n  1 

i+1 

if 

3  =  J+i. 

o  >  o, 
n  1 

0 

if 

t  <  tj 

,  n  *  a 

It] 

if 

t  <  tj  , 

n  <  a- 

(“n^R*  if  tj  <  t  <  tJ+1  , 


y  „[t]  = 
*R 
n 


0  if  tj  <  tj 


<uR>un+1*  fcJ  <  <  fcJ+1 


L<n  i  {»•  8  0  <  k  <  n} 


l  =  u*;,**}  . 

n  n  n 


:  n  <  k  <  a-1> 


Define 

(4.13) 


L  =  U  e  L<nu  Ln  :  Tt(tJ+)  f  0} 


For  l  e  L  let  o'?,  q  »  L,  R,  be  defined  by 


and  aet 
(4.14) 


ui(tj+)  “  T(a1l  V  • 


5  U  eL‘  !  °M*  (°X>} 


*  L  R 

Moreover  if  l  G  L  and  G  (0^,0^),  then  let 

.  v4 

0  =  TT  ' 

V°* 


L*  a  «■  »  *  •  **'e  *  .  ■  ,  «  .  •  . 


•».vV.-:.>- ,v_\ 
■/  V  vv  .- 


!-  ■•-  /..v.-  .'  .-  V  .1  .■  ■  ■  v  A-\*  •.«V  .V\. 


and  set  (cf.  (4.8)  -  (4.10)) 

(4.15)  :  {*’  i  1  6  L*,  0M  e  (<^/0^).  q  «  L  or  R)  . 

Finally,  we  define 

(4.16)  W^(j+i)  =L^m u  LMuL>nuL*  , 
in  the  case  Y?_  f  0. 

1J 

Since  y ^  is  defined  for  all  l  e  ( J+ 1 )  this  completes  the  definition  of  Mp(J+1> 

and  r^(J+1) .  we  leave  it  to  the  reader  to  verify  from  the  above  construction  that 
Properties  (4.1)  -  (4.3)  are  satisfied  by  Mp(J+1)  and  r^(J+1).  This  completes  the 
definition  of  the  approximate  characteristics. 

We  say  that  Yt  e  T(J)  is  cancelled  at  (x^tj)  if  i(j)  «  i  and  j  -  )£• 

(4.16)  it  is  clear  that  Yj  e  T(J+1)  is  cancelled  at  (x^,tj)  if  and  only  if 
i  e  L>  u  L,  or  i  -  iR  where  n  ■  q. 


From 


§5.  PROPERTIES  OF  APPROXIMATE  CHARACTERISTICS 

In  this  section  we  study  properties  of  the  characteristics  T(J)  and  index  set  M (J) 

defined  in  the  previous  section.  Let  uh  be  a  given  approximate  solution  generated  by  the 

random  choice  method  from  initial  data  uq.  The  sets  T(J)  and  M  (J)  associated  with 
uh  are  determined  by  the  choices  of  uQ,  J,  h  and  T,  which  we  take  to  be  given.  We 
let  M  =  M(J),  P  =  T(  J) .  Throughout  the  remainder  of  the  paper,  is  taken  to  swan  £ 

L  tel 

whenever  L  £.  M-  We  now  develop  some  definitions. 

Let  N(t)  denote  the  index  set  for  the  "null"  characteristics  that  are  either 
cancelled  or  else  are  created  by  "nonlinearities*  in  time  [0,t],  t  <  tj*  i.e., 

(5.1)  N(t)  =  U  6  M  :  >  0  or  t’  <  t}  . 

We  partition  N(t)  into  N0(t)  and  N^(t)  as  follows: 

(5.2)  N0(t)  =  U  e  N(t )  :  t°  -  0}  , 

(5.3)  N^(t)  S  {i  e  N(t)  :  t°  >  0}  . 

By  Property  (4.3),  the  index  set  M  partitions  into 

M  «  M*  u  M 

where 

M  +  s  U  e  M  :  sgn(i)  >  0}  , 

=  (tell  :  sgnlt)  <  0}  . 

We  call  M+  [reap  M-1  the  set  of  rarefaction  wave  [reap,  shock  wave]  characteristics. 

For  p  =  1,...,n  and  q  m  +  or  ~,  define 

(5.4)  Mq  i  Mp  n  Wq  • 

Similarly,  define 

Nq(t>  =  N(t)  n  Mq  , 

and  set 

Nqn(t)  =  Wq(t)  n  Nn(t)  , 
po  p  u 

=  Nq<t)  n  Hu(t)  , 
pp  P  P 

so  that  N^(t)  is  the  disjoint  union  of  NqQ(t)  u  N^tt). 
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We  say  that  two  characteristics  Yjj  and  Ym  intersect  at  time  tj  if  l(j)  “  m(j) 
and  i(j-1)  j*  m(j-1).  The  following  lemma  which  can  be  verified  by  induction  on  J 
implies  the  uniqueness  of  intersection  times  for  characteristics  in  different  families. 


LEMMA  (5.1).  If  t  e  Mp  #  m  e  Mp  ,  and  p  <  p* ,  then  xm  -  is  a  nondecreasing 

function  of  time  for  t  e  [t®,th  n  [t^.t1]. 

k  x.  mm 

The  next  lemma  follows  from  Properties  (4.1)  through  (4.3): 


LEMMA  (5.2).  If  i,  m  e  Mp  ,  sgn(f)  /  sgn(m),  and  i  intersects  m  at  time 
tj  <  t  <  tj,  then  either  l  e  N(t)  or  me  N(t).  Moreover,  if  I  and  m  are  both  shock 
wave  characteristics,  then  there  is  at  most  one  intersection  time  tQ,  and  for  all 
t  e  (tQ,tJ)  n  [t^t^)  we  have 


Define 


Ni;ja)  = 


xI<t>  “  *m(t) 


|Yt<tj+,|  -  |Yi(tj-)|  if  Mj)  -  i  , 


otherwise 


We  call  N^jft)  the  nonlinearity  contributed  by  the  characteristic  y  ^  at  the  mesh  point 
(x^.tj)  (cf.  16  ]•  In  particular,  since  Yt(t)  18  constant  on  (t£,t j] ,  we  must  have 
Nj^jU)  »  0  unless  4  e  N(tj)  and  either  j°  -  j  or  -  j.  Define 


c^U) 


|TTJt(tj-)|  if  4(j>  -  i  and  j  -  jj 


otherwise 


IYjjU.-*-)!  if  I(j)  -  i  and  j  -  j° 


0  otherwise  , 


so  that 


Nij(i)  -  -Cij(4)  +  Ei3(£)  . 


We  call  Ci3(t)  the  cancellation  and  E^tt)  the  error  in  the  nonlinearity  N^U).  The 
following  proposition  is  a  restatement  of  Lemma  (GL1)  in  the  language  of  approximate 
characteristics : 


-  **.  -*.  .*.  ••  v.  .*. 

'j  *  ’  *  •  ,  1 


v>v-v- v>v-v. 


v  %• >,■  . 

•  »  m  ■  fc  m  * 


\.Y  A  A  A  a  A*  .n' 


PROPOSITION  (5.3).  There  exists  GQ  >  0  such  that 


(5.10) 

and 

(5.11) 


l 'd''1  ‘ 


)l£,  v«’l-  ^1  “vi, 


for  p  *  1,...,n  and  q  «  +  or  -. 

(Note  that  since  E^lt)  “  C^_ ^ ( E. )  -  0  for  t  fL  N(tj),  sums  over  M  and  can  be 

replaced  by  sums  over  N(tj)  and  A/p(tj)  in  (5.10)  and  (5.11),  respectively.) 
COROLLARY  (5. 3).  There  exists  Gq  >  0  such  that 


(5.12A) 


(5. 12B) 


l  |Ytl  <  G0[0(0)  -  Q(t)l 


yti 


1  htl  <  Govo 


Proof*  If  t  6  then  tj  >  0*  Thus  by  (5.5), 


htl  -  ■  o  oU) 


and  so  by  Lemma  (GL1 ) 


V*’ 


l  Iy.I  <  l  l  E  .  (Jl>  <  l  G  D  <  G0[Q(0)  -  Q(t)) 

iJM  3  ij  3 


This  verifies  (5.12A).  For  (5.12B),  note  that  since  |  Yj^(  t )  |  is  constant  on  [t£,tj|), 
must  have 

1  Ini  <  Govo  • 


Thus  by  (5.12A) 


l  |Tt! 

M 


I  htl  +  1  |Ytl  <  G0V0  . 


M'%  <V  Mf(tJ) 


For  j  <  J,  define  the  set  of  indice  pairs  associated  with  waves  that  approach  at 
time  t  (cf .  (3. 11)1 i 

Aft)  -  A’(t)  u  A2(t) 


where 
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1 


(5.13) 


A1  < t )  =  {<4,m>  :  t,  m  6  Mp  for  some  p  -  1,...,n,  x^tt)  <  x^Ct)  , 
and  either  4  or  m  Is  in  M*}  , 


(5.14)  A2(t)  =  {<4,m>  s  4  e  M  ,  a  e  U  ,,  for  p  <  p',  and  x.(t)  >  x  (t) } 

P  p  k  in 

(Here  we  use  <l,m>  to  denote  the  set  or  unordered  pair  Note  also  that 

X4(t)  <  Xm(t)  implies  xA(t)  +  t  /  x^t) ,  so  that  |Y4<t)  |  ?  0  +  | Ym(t )  | ) . 

1  2 

He  call  A  (t)  and  A  (t)  the  index  seta  for  the  characteristics  that  approach  at 
time  t  and  are  in  the  same  and  different  families,  respectively.  Define 

(5.15)  A^j  =  {<4,m>  :  4(j)  -  m(j)  -  i  and  <4,m>  e  A(t^-)} 

The  following  is  a  consequence  of  Property  (4.3)  together  with  the  definition  of  Q  and 
Dij: 


LEMMA  (5.4): 

He  have 

(5.16) 

Q(t)  -  I  1**1  | Ym 1  . 

<4,m>8A(t) 

and 

(5.17) 

Dij  “  I  Iy4I  1 Yml  • 

J  <4,m>eAlj 

From  here  on 

out,  we  write  £ 

in  place  of  l  ,  etc. 

A 

<4,m>eA 

Define 

A0(t)  and  A^(t) 

as  follows: 

(5.18) 

A0(t)  =  {<4,m> 

e  A(0)  :  either  4  or  m  is  in 

«♦ 

*-r-> 

(5.19) 

A|,(t)  5  {<4,m> 

e  A(t)  :  either  4  or  m  is  in 

yt» 

LEMMA  (5.5):  Let  A'  c  A(0)\A(T),  T  <  tj.  Then 

(5.20)  A(T)  =  A(0)  U  A|(t)\A'  • 

Proof:  By  the  definition  of  A(t),  <4,m>  e  A(t)  implies  that  Y4(t)  t  0  and  Ym(t)  +  0, 
thus  A(T)  n  A'  “  0.  Assume  then  that  4  and  m  are  in  M\Njj ( T)  and  <4,m>  e  A(T).  It 
remains  only  to  show  that  <4,m>  e  A(0).  But  by  the  definition  of  A,  the  only  way  <4 ,m> 
could  fail  to  be  in  A(0)  Is  if  4  and  m  intersect  at  two  distinct  times  in  [0,T], 
which  contradicts  Lemma  (5.1)  or  (5.2).  For  example,  if  4  e  Mp  and  m  eMpi,  p  <  p' , 
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(5.26) 


fi(0)  -  Q(T)  >  L2 
G0 

Since  T  is  fixed  here,  we  set  A/*1  =  hft(T),  etc. 

P  P 


Proof j  Assume  that  (5.25)  holds.  Then  for  some  p  e  {l,...,n}  and  q  «  +  or 
must  have 
(5.27) 


I  1**1  ‘ 

N?(t) 


For  convenience,  assume  q  «  +.  By  (5.7)  and  (5.8), 


(5.28) 


i  1**1  «  l  l  *•**<*>  +  L  l  C  it) 

»;,«»  ‘i  iAti  1  « 1 

P  P  p 

<  l  GoDij  +  l  l 

ii  J  il  +  J 

Np(t) 


where  the  sum  on  i,  j  is  over  {(i,j)  :  -»  <  i  <  t.  <  t}.  But  V  d  < 

j  ij  ^ 

”  Q<t)],  together  with  (5.27)  and  (5.29),  implies  that  for  some  G0  >  0, 

i  5.29)  l  l  cij  (*.)  >  ^  -  G0(Q(0)  -  Q(t)] 

13 


and  so  by  ( 5. t 1 ) 

(5.30) 

Thus  we  can  apply  (5.11)  again  with  q  «  -  to  obtain 


l  °ij  >b-  8o«,0>  • 


I.  I  C^U)  -  I  _  <  G0  [Q(  0)  -  Q(t)]  , 


ij 


N  <t) 
p 


ij 


'ij 


or 

(5.31) 


l  l  CijU)  >  IS  '  «30tQ(0)  -  Q(t)l  . 


W*(t, 


Therefore,  by  (5.7)  and  (5.8), 


l  \ye\  >  l  1  c^ts)  - 1  i  E^ts) 


sewp(t) 


seAl  (t) 

P 


ij 


sew  (t) 

p 


ij' 


>  %  -  vc(0)  ■  s(t)J  • 


we 


-32- 


He  now  have  that  for  some  Gq  >  0,  both 


C5.32A) 


l  I  Yr  |  >  -  G0[Q(0)  -  Q(t)] 


(5.33B) 


l  lYs!  >  "  Go > ]  . 

aeH+ 

p 


By  Corollary  (5.3),  this  implies  that  for  some  GQ  >  0,  q  =  +  and  -,  we  have 


(5.34) 


l  l*tl  >  ^  -  GqLQ(O)  -  Q(  t ) )  . 


Let  A'  5  {<r,s>  :  r  e  M^q,  s  e  N^q).  Then  by  (5.2),  (5.5)  and  the  definition  of  A, 


A'  A(0)\A(T).  Moreover, 


I, \yt\  |Tj  >  l  l  |frl  !y9I  >{^“G0[Q(0)  -Q(t))}2 
A '  reA/+.  sew" 

p0  pO 


>  -=-j  -  G0(Q(0)  -  Q(  t )  ] 
4n 


where  we  use  the  fact  that  L  <  VQ  <  GQ  and  tQ(0)  -  Q( t ) )  <  GQ  for  some  G0  >  0.  Thus 


by  Lemma  ( 5.7) , 


Q(0)  -  Q(T)  >  — '—j  L2  -  [Q(0)  -  Q(T>]  , 

GQ4n 


(5.35) 


Q(0)  -  B(T)  >  - '—j  L2  . 

8G0n 


Since  Gq  is  generic  we  can  take  - -  to  be  Gq,  so  (5.35)  establishes  Proposition 

8G0n2 

(5.8). 


Corollary  (5.8).  There  exists  G0  positive  so  that  if  T  <  tj  and 


(5.36) 


0(0)  -  Q(T)  <  L  , 


I  |Yj!  <  g0l/2  . 
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We  next  define  index  sets  for  characteristics  that  intersect  a  given  characteristic  in 
a  given  time  interval.  Let 


(5-38) 


(5.39) 


8  <  j)  =  (m  e  M  :  m<  j)  *  £(  j)  and  <4,m>  e  A^(t  .-)}  , 

*  3 

Bi<5>  s  _  l  K\  ' 


and  for  t1  e  [t^ tj  1 )  and  t2  e  ttj2'tj2+1)'  let 

j2 

(5.40)  B4[t1#t2l  =  l  Bjlj) 

jl 

In  addition,  for  i  e  M“,  define 


®  AR«  3  J  H  {meWp  !  m{3)  "  *<3)r  m(j-1)  >  i(j-1)} 


8^(3)  5  {meM~  :  m(j)  -  t(j),  m(j-1)  >  i(j-1)}  , 

B-,(3)  -8;r(3)u  B-l(3)  • 


For  q  •  L,  R  or  absent,  define  8Jq[t1,t2]  as  in  (5.40).  (Note  here  that  script 
denotes  a  set,  while  upper  case  B  denotes  a  real  number.) 

In  the  next  section,  we  use  the  following  technical  result: 


PROPOSITION  (5.9).  Let 


(5.41) 


w:uL,iRi  =  «®>  e  iw} 


Assume  that  s  e  M"\N(t),  s(0)  =  iL,  that 


(5.42) 


and  that 


(5.43) 


l  UJ  >  4L 

VW 


Y*  <  L 


for  all  i  e  tiL,iR)  such  that  y|q  <  0.  Then  there  exists  G0  >  0  such  that,  if  L  <  u 


(5.44) 


1  IyJ  >  L  , 


BpR[0,T] 


(5.45) 


Q(  0 )  -  Q(T)  >  G"'LZ 
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Moreover,  if  the  hypotheses  are  satisfied  with  s(0)  =  iR  and  8“l[o,TJ  >  L,  then  again 
(5.4S)  follows. 

To  prove  Proposition  (5.9),  we  use  the  following  lemma: 


LEMMA  (5.9).  Let  {L^}®^  be  a  nonnegative  sequence  of  numbers  which  satisfies 


(5.46) 


0  <  L  <  L  , 


(5.47) 


l  L.  >  3L  . 
i=1  1 


(5.48) 


V  L.L.  >  L2  • 
1<i< j<m  1 


Proof .  Let 


Therefore, 


i*  =  inf  ti  :  V  L.  >  l} .  Then  L  *  <  L  implies  that 
i=1  1  i 

l  \  >  M  . 


i«i  +1 


i  LL  >  (  I  L )(  L  )  >  L2 

<j<m  J  i»1  i*i  +1 


Ki<j<m 


Proof  of  Proposition  (5.9):  Let 


8o  5  B“r[0,T]\N(T) 


8.  i  8  ro,T]  n  N(t)  . 
P  pR 


L0  *  I  ' 


*  l  1^*1  • 

'  B„ 


Without  loss  of  generality,  assume  Lg  >  3L.  If  not,  then  L^  >  L,  in  which  case 


Proposition  (5.8)  implies  that 


C(0)  -  Q(T)  >  L2  , 
G0 


* •* •-*'  y 
_  .*  -* 


-•.y-j 
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which  gives  (5.45).  Let 


and  define 


8io  2  U  e  80  :  i(0)  =  i}  , 


Li  “  I  I’fjjl 

B  * 
iO 

Let  i"  *  max{i  :  L^  j*  0},  and  let  i'  «  min{i",iR).  By  (5.42)  and  (5.44), 


1  L.  >  3L  . 


i«i 


Moreover  by  Property  (4.3),  L^  “  yfg  <  0  and  L^  «  0  otherwise,  so  that 

L^  <  L  for  all  i  6  [iL,i*].  Thus  the  conditions  of  Lemma  (5.9)  hold,  and  we  have 

l  L  L  >  L2  . 

i  <i< j<i'  J 
L 

But  by  definition,  <t,m>  e  A'  (0)\  A'(T)  if  l,  m  e  8g  and  4(0)  ^  m(0).  Thus 

A(0)\  A(T)  >  l  It, I  I Y  |  >  l  L.L,  >  L2  • 

A*  i  <i<j«i'  1  3 


Therefore,  by  Lenina  (5.7)  we  conclude 


Q(0)  -  Q(T)  >  7T-  L2 
^0 


Finally,  define  X^lt] ,  the  speed  of  the  characteristic  Yt  at  time  t  as  follows: 
Let  4  e  Mp  ,  t  e  [tj,tj+1)  n  (t“,t3),  and  assume  that  4(j)  -  i.  Then 


X4[t]  2 


X(  Y*\  )  if  t<0  , 


X  (u,  [tj  )  if  l  >  0 
P  * 


We  now  estimate  the  change  in  u£,  u£  and  X^  for  i  eM(J)\N(J) 


LEMMA  (5.10).  Let  4eM(J)\N(J).  Then 
(5.50) 


u^ttj+]  -  »2ltj-]|  <  G0<Bi(j>  ♦  Dl(j)/j} 


(5.51)  |X4(tj+]  -  X4[t.-]|  <  G0{Bi(j)  ♦  B;(j)  -  Ct(j)fj  + 

for  q  «  L,  R. 
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COROLLARY  (5.10).  Let  4  S  M\N(T),  t-j  <  t2  <  T  <  tj.  Then  for  q  =  L  and  R, 

(5.52)  lultt21  <  G0{Bt[t1't21  +  _  ' 

(5.53)  lXilt2]  "  Xftt1]l  4  G0tB4It1*t21  +  Bi[t1,t21  +  lClt1>  "  Q<t2))1/2}  • 

Moreover, 

(5.54)  Iu?tt2l  -  ujtt,]  <  G0Vq 

(5.55)  1 XA  [t2]  “  W  I  <  G0V0  ‘ 

Let  a  be  a  fixed  equidistributed  sequence,  0  <  a^  <  1  115,  16]  .  Let  I  = 

(c,d)C[0,1]  and  let  0  <  <  ^2*  Define 

N(X,N1,N2)  =  Card{ j  e  [N^Nj]  :  a.  e  i)  . 

The  following  is  a  result  regarding  equidistributed  sequences  [15,  16] s 

LEMMA  (5.11).  For  every  M  >  1  there  exists  N0  >  1  such  that,  if  N  >  N0,  then 


for  all  j  <  MN,  I  =  (c,d)C[0, 1] ,  |l|  =  |d-c|.  Moreover,  for  the  best  equidistributed 
sequences , 

(5.57)  N(M,a>  <  GqM'1  , 

where  N(M,£)  denotes  the  infemum  of  all  such  N0  for  a  given  a. 

For  4  e  M,  t  e  ItG,t^),  define 

(5.58)  E^(t)  *  x4(t)  -  x^ft®)  -  JtQ  X^ItJdt  . 

THEOREM  (5.12).  For  each  equidistributed  sequence  a  there  exists  a  positive  function 
lj(U)  <  1,  j(()J)  ♦  0  as  u  ♦  0,  such  that  the  following  two  statements  hold: 

If 

(5.59)  h  <  |{(U)v,  t  e  [t°,t’>, 
then 

|y4l  <  V 

for  all  4  e  M+(J)  such  that  t^  =  0,  t^  >  v. 
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Since  z^(u^[0]  )  -  z^fu^tO])  is  positive  for  i  e  M+  and  negative  for  l  e  M~;  and 

1-6  — 

since  z increases  from  m  to  —  going  from  x^  to  x  at  time  t  =  0,  we  can  write 

<»■»>  Y', • 

Since  Rf.Ro  u  U  n  N(T)},  (6.57),  (6.58)  and  (6.59)  imply 

(6.60,  RU0  h  3  t0' 


where  X  is  a  set  of  small  measure. 


(6.61) 


m{x}  <  iT  <  1SR  ' 


where  we  have  used  (6.22)  to  estimate  6. 


By  (5.59),  if  re  Rg,  then 


(6.62) 


IV  <  16M 


because  h  <  G„^(tss)^  by  (6.14).  Thus  (6.60),  (6.61)  and  (6.62)  immediately  imply  the 

U  1  oM 

existence  of  r  e  R0  such  that  (6.46)  holds. 

Similarly,  we  conclude  as  in  (6.60)  that 

UI1  -  l0‘  ii,Nx  ' 

3o 

where  X  satisfies  (6.61),  which  directly  implies  the  existence  of  s  e  SQ  satisfying 
either  (6.47),  (6.48)  or  (6.49),  (6.50). 

Define  [cf.  (5.38)  -  (5.40)] 

(6.63)  Rjj  =  U  e  R0  t  x4<0)  <  xr(0)}  , 


(6.64) 


S*  2  u  e  S0  :  x4(0)  >  xs(0)>  , 


=  {l  e  Mpt  :  p1  /  p  and  intersects  yg  in  [0,  T/2]} 


(6.67) 


B0  =  i  Utl 

B0 

8r  s  8r[0,  T/2],  Br  2  l  ( y £ 


(6.69, 


6S  S  8g[0,  t/2],  bs  s  l  |y4| 
B 

r 


[0,  T/2],  B"  i  l  jyj,  q  =  L,  R  or  absent. 


>V-W 


• '  •  *  *v>  • »  *  »'• , 


ind  either 

6.47) 
ind 

6.48) 

>r  else 

6.49) 


VUsl°])  4  In 


|yp  I  >  - 

1  i  O1  M 

8 


Zp(us 1  e  *8M  '  4M^ 


ind  there  exists  iL,  iR  with  iA  <  iL  <  ifl  <  iR  <  iB  such  that 


6.50) 

or  all  i  eMp[iL,iR],  and  [cf.  (5.41)) 


s 

l^ol  «l 


6.51) 


6.52) 


I 


M  [i.  ,i  ] 

p  L 


1  ltj>^ 


Mp[is'V 

"roof .  We  first  verify  the  existence  of  [  e  Kg  satisfying  (6.46).  Let 
6.53)  L  =  U  e  M  i  x4(0)  e  [xA,x]}  , 

P  *'  J  Zp( 


,nd  let  1^  denote  the  interval  [z  (u£[0]  ),  z  (u^[0))).  Let  1 1^  [  denote  the  length  of 


he  interval  1^.  Let 


m  5  Inf {z(u  (x) )  s  x  e  (x  ,x]  } 
0  A 


6.54) 

ly  (6.39)  and  (6.41), 

6.55) 

Toperty  (4.3)  states  that  the  approximate  characteristics  partition  the  waves  in  the 
ipproximate  solution  uh  at  each  time  step.  This  together  with  (6.40)  implies 

6.56) 

iy  Leiana  (6.4), 

6.57) 


-M  <  m  <  |  . 


I  ntl  . 

L\  Un 


nd  since  we  assume  Q(0)  -  y(T)  4  - Corollary  (5.8)  implies 

<GM)2  . 

6-58)  l  |l  |  <  |  . 

LnN(T)  1  M 
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s  *. 


the  case  when  5*  *  -Inf  {2  (u^J(x))}  being  handled  analogously. 

M  D  U 


By  (6.39)  there  exists  a  point  x  =  (i  +  -j)h  such  that 


(6.40) 

Define  iA  and  iB  by 


Zp(U0(x)) 


ift  =  Sup{i  :  i  <  i  +  ~  and  z  (u^fx^))  < 


(6.41) 


iB  =  Inf (i  :  i  >  i  +  j  and  *  <Vxi  +  ,)  <  ^ 


Let 


XA  2  UA  "  2)h  ' 


(6.42) 


XB  “  UB  +  2)h 


Define  R  M”  and  S  c_  M*  by 


R  i  {i  eW  :  x.(0)  e  (x.,x)>  , 


(6.43) 


S  =  {t  e  M  i  x  (0)  e  (x,x  )}  , 

p  X  2 


and  let  [cf.  (5.1)] 

(6.44) 


RQ  =  R\W(  T)  , 


SQ  =  S  \W(T) 


LEMMA  (6.5).  The  following  estimate  holds: 

(6.45)  |xB  -  x.|  < 


,  -  eM 

XB  -  xAl  '  T  • 


Proof.  By  (6.40), 


E  *  /UU£(x)ldx  >  &  lxB  -  XA 


Solving  for  |xB  -  xA|  gives  (6.45). 


PROPOSITION  (6.6).  If  Q(0)  -  Q(T)  <  — — j,  then  there  exists  r  e  R0  and  s  e  S0, 


(GM) 


s(0)  «  is,  such  that  the  following  conditions  hold: 
(6.46) 


Up(«J(0J).  *p(u*[0]))c  (if.  2_,  , 
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/  s 


He  use  this  to  construct  two  characteristics  Y  e  T  and  Y_  6  r  which  would  intersect 

r  p  s  p 

before  time  T  if  there  were  no  interactions  to  deflect  the  speeds  of  these 

characteristics  from  their  initial  speeds  at  time  t  =  0.  We  then  assume  that  Q(0)  - 

Q(T)  <  — r.  By  this  assumption,  there  exists  such  Y_»  Y_  which  are  not  cancelled  in 
<GM)2  r  8 

time  [ 0, T] •  We  then  use  the  same  assumption  to  obtain  estimates  for  the  change  in  the 

speeds  of  Yr  and  Ys  between  time  t  =  0  and  t  =  T.  These  estimates  are  sufficient  to 

guarantee  that  in  time  [0,T]  the  wave  speeds  of  Yr  and  Yg  are  not  deflected  enough 

from  their  initial  speeds  to  prevent  them  from  intersecting  before  time  T.  This 

intersection  implies  that  one  of  them  is  cancelled  before  time  T  [cf  Property  (4.1)]. 


This  is  a  contradiction  and  thus  we  conclude  that  Q(0)  -  £>(T)  > 


We  first  use 


Lemma  (6.1)  together  with  (6.30)  or  (6.24)  to  estimate  the  total  variation  in  zf 
contributed  by  characteristics  not  in  the  p-family. 


LEMMA  (6.4).  Our  assumptions  on  u"  imply 


(6.36) 


1  |.  (uJiod  -z  (ujionl 

M\Mp  p  p 


Proof :  A  restatement  of  Lemma  (6.1)  in  the  language  of  approximate  characteristics  is  that 


(6.37A) 


|zp(u*[0])  -  Zp  ( u  ^  [  0  ]  )  |  <  I  Ya(0)  I 


for  all  A  e  M\Mp;  and  if  z  is  a  coordinate  system  of  Riemann  invariants,  then 
F  G 

(6.  37b)  Izp(u>l)  -zp(u>))|  <-||Yt(0)|  . 

M 

But  by  Property  (4.2) 


(6.38) 


l  |t.(0)|  =  v  4  V 


Thus  statement  (6.36)  follows  directly  from  either  (6.37A)  or  (6.37B)  by  estimating  the 
right  hand  sides  using  either  (6.30)  or  (6.34A),  respectively. 

By  (6.26)  there  exists  p,  1  <  p  <  n,  such  that 

Sf  -  Sup  |z  <u£(x>)|  . 

X 

For  convenience  we  assume  that 


(6.39) 


4  =  Sup  {z  (u*?(x) )}  , 

M  p  0 

X  r 


I  •  •  •  -  «  «  •  *  *  *  ,  •  '  *  *  •  *,  •  •  ,  * 

V  V  */  ■>  •>.*  V  V  *  *  v  •  * 


•  V  *  .  '»  *  »  '  *  >  *  *  • 

■-w - 


-  V 


The  choice  of  zp  for  the  p-wave  strength  parameter  determines  the  definition  of  wave 
strength  for  the  characteristics  in  T  [cf  (4.6)].  In  this  case,  for  t  e  Mp, 

5  Vufltl)  “  zp(uict3) 

defines  the  signed  strength  of  the  characteristic  y^  at  time  t  <  tj.  Recall  that 
Yt<t)  =  Yt  is  constant  for  t  e  [t°,t^),  and  identically  zero  elsewhere. 

The  idea  in  the  proof  to  follows  is  this:  since  lu0(*)lg  *  ~  and  I uQ < * ) * 0  “ 
there  must  be  a  "spike”  in  the  initial  data  of  height  -j  and  width  on  the  order  of  eM. 


In  the  analysis  to  follow,  it  is  important  to  bound  changes  in  Xp  by  changes  in  Zp. 

This  can  be  done  because  z  ♦  u  is  a  regular  map,  and  71  *R  >  0.  Define 

P  P 


31  (T  (o»u  ) ) 


(6.20) 


1  =  Inf 


where  a  =  z^(u)  -  Zp(uL)  and  the  infemum  is  taken  over  all  values  of  o,  uL  such  that 

u  =  T  (cj;uL)  6  U. 

P 

Now  let  Gq  be  large  enough  to  satisfy  all  previous  conditions,  as  well  as 


(6.21) 


Gq  >  Max{l,  l”  ,  V)  . 


Define  the  following  constants: 


(6.22) 


«  =  «<G  )  =  , 

64G„ 

0 


(6.23) 


G  =  G(6,G.)  -  , 

0  6"* 


(6.24)  M2  =  M  (6,Gg)  “  -j—  • 

Choose  U  to  be  a  sufficiently  small  neighborhood  of  u  “  0  so  that  (6.14)  holds  and 


(6.25) 


U  £  {u  :  |  u  |  <  M^ 1 } 


for  i  «■  1,2.  We  prove  the  follow  theorem  which  is  a  restatement  of  Lemma  (2)  of  the 
Introduction. 

THEOREM  (6.3):  Let  uh  be  a  fixed  random  choice  method  approximate  solution  satisfying 


( 1 ) .  Define 


ib  ,uo<°>'s  ' 


(6.27) 


e  =  lu“(0)l  ,  , 

L 


(6.28) 


T  =  E(GM)  . 


In  the  case  that  z  is  not  a  coordinate  system  of  Riemann  invariants,  assume  that 


(6.29) 


Vn  <  V  , 


where  V  is  sufficient  for  U  in  Lemma  (GL2)  and 


(6.30) 


’‘3  • 
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(6.13) 


(*p.  CuH>  -  zp.(uL,|  <  ®  |^|  . 


Proof ;  Statements  (6.12)  and  (6.13)  follow  directly  from  (6.3)  together  with  the  fact  that 

Sp(uL)  makes  P.W.  contact  with  Rp(uL)  at  the  state  uL  (c f.  (2.1)]. 

It  is  clear  from  (13)  that  if  U  is  sufficiently  small,  then 


(6.14) 


|zp(uR)  -  Zp(uL) I  >  2|zp.(uR>  -  zp.(uL)| 


for  all  p-waves  with  left  and  right  states  uL  and  uR  in  U,  p'  ^  p. 

The  next  lemma  is  a  technical  but  elementary  uniform  estimate  for  the  speed  of  a  p- 

shock  in  terms  of  zp. 

LEMMA  (6.2):  Let  S  denote  a  p-shock  with  speed  s  and  left  and  right  states  uL  and 

D 

u  .  Then  there  exists  a  constant  M-)  >  0  and  a  constant  Gq  >  0  depending  only  on  M1 
and  f,  such  that,  if  M  >  M.  and 


(6.15) 


|«q|  <  £  ,  q  *=  L,R  , 


(6.  16) 


z  (uL)  +  z  (uR) 


16M  ' 


(6.17) 


for  all  p*  f  p,  then 


lv<u)  ■  V(u  >!  ‘  V  ' 


(6.18) 


X  (u)  -  a  >  — 
p  G  M 


Proof :  Lemma  (6.2)  expresses  in  a  uniform  way  the  fact  of  Lemma  (2.1)  that 


(6.19) 


*  0(1)  |B|* 


The  existence  of  follows  from  (6.3)  and  the  uniformity  follows  from  the  compactness 

of  U  together  with  the  fact  that 


~  X  >o 
9z  p 
P 


satisfying  (6.1)  together  with 

(6.5)  VQ  =  1  |y*0j  <  V  . 

i/P 

takes  values  in  U  for  all  x  and  t,  and  satisfies  the  conclusions  (3.18)  -  (3.20)  of 
Leona  (GL2 )  as  well.  He  say  that  M  is  sufficient  for  U  and  V  in  Lemma  (GL3)  if 

(6.6)  V0  <  V  , 
together  with 

(6.7)  'ujj( ‘Jig  *  ”  . 

M 

guarantee  that  u*1  takes  values  in  U  for  all  x,  t,  and  satisfies  the  conclusions 

(3.18)  -  (3.23)  of  Lemma  (GL3)  as  well.  It  is  clear  from  the  statements  of  Leonas  (GL2) 

and  (GL3)  that  for  every  neighborhood  U  of  u  •  0  there  is  a  V*  such  that  any 

V  <  V'  will  be  sufficient  for  U  in  Lemma  (GL2)>  and  in  the  case  of  a  coordinate  system 

of  Riemann  invariants,  for  every  V  and  V  there  is  an  M'  such  that  any  H  >  M'  will 

be  sufficient  for  U  and  V  in  Lemma  (GL3 ) . 

For  u  @  U,  let  ju|  =  Sup|z  (u)|,  and  for  functions  u  :  R  ♦  0,  define 

P  P 

(6.8)  |u(x)|  =  Sup|s  (u(x))|  , 

P.x  P 

(6.9)  lu(*)l  .  a  J^Jutx)  |dx  , 

L 

(6.10)  lu(-)lc  =  Sup{ | u( x) | }  . 

s 

X 

The  following  two  leonas  will  be  needed.  The  first  lemma  estimates  the  change  in  zp, 
across  a  p-wave,  p'  /  p: 

LEMMA  (6.1):  Let  be  any  p-wave  with  left  and  right  states  uL  and  uR  satisfying 

(6.11)  |u*|  <  ^  ,  q  -  L,  R  . 

Then 

Q 

(6.12)  |zp.(uR)  -  zp,(uL)|  <^1^1  • 

Moreover,  if  z  is  a  co  rdinate  system  of  Riemann  invariants,  then 
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$6*  THE  MAIM  ESTIMATE 

In  this  section  we  study  approximate  solutions  uh(x,t)  generated  by  the  random 


choice  method  from  initial  data  Ug(x)  which  satisfies 

(6. 1)  u£<±  »)  =  0  . 

We  study  the  approximate  solutions  in  a  coordinate  system  of  Riemann  invariants  if  one 
exists,  and  if  not,  then  in  a  coordinate  system  that  is  a  good  approximation  to  a 
coordinate  system  of  Riemann  invariants  near  u  -  0. 

Thus,  let  z  =  (z1,...,zn)  denote  a  coordinate  system  of  Riemann  invariants  if  one 
exists;  i.e.,  in  this  case  assume  that  the  mapping  u  z  is  a  1-1  smooth  map  taking 
0  +  0,  and  which  satisfies  the  condition 

<6*2>  Tz-  * 

k 

Such  a  coordinate  system  exists  if  and  only  if  there  is  a  choice  of  eigenvector  fields 


such  that 


[Rj,R£]  -  0 


for  all  j,  k  e  {l,...,n},  where  [  ]  denotes  the  Lie  Bracket.  A  coordinate  system  of 
Riemann  invariants  always  exists  in  the  case  n  “  2. 

If  a  coordinate  system  of  Riemann  invariants  does  not  exist,  then  choose  z  « 


( z1,...,zn )  to  satisfy  [cf.  6,  17] 


IT  *  Rk(0) 


In  either  case,  let  zp  be  the  wave  strength  parameterization  of  the  p-shock-rarefaction 
curve  in  a  neighborhood  of  u  *  0  [cf.  (2.1)]. 

For  example ,  we  take  o  =  zp(u)  -  zp(uL)  in  the  equation  u  *  T^ta;^),  so  that 
(6.4)  yP  =  zp(uR)  -  zp(uL) 

defines  the  signed  strength  of  a  p-wave  1P  with  left  and  right  states  uL  and  uR. 

We  let  U  denote  a  neighborhood  of  u  -  0  in  which  Riemann  problems  are  uniquely 
solvable  such  that  zp  is  a  regular  wave  strength  parameter  for  all  p-wave  curves  in  U, 
and  such  that  Lemma  (GL1)  (cf.  (3.15),  (3.16)]  applies  with  this  measure  of  wave  strength. 
We  say  that  V  is  sufficient  for  U  in  Lemma  (GL2)  if  any  approximate  solution  uh 


■>vv>:v 


•  .**  ■-  \  >.  •_  - 
•  #  a  t  •  •  ,* 

r  •.*  •  ’  %  o  -  ■  v  '* 


mm  *«  N*  *,*  *  * 

-  **.  *’s  *\  .*.  .“.  **  . 
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.*»**.*■ 


thus  verifying  (5.60).  For  the  best  equidistributed  sequences,  N(M,a)  “  GQy  1  by  Lemma 
(5.11),  and  so  in  this  case 

(5(0)  -  g/ 

by  (5.63)  and  (5.64).  This  completes  the  proof  of  Theorem  (5.12). 


V*. 


N 

|Et(1,N)|  -  |  l  E4(j)| 

j=1 

<5.71)  <  <£  X^hN,  -  <1  -  f  xl)hir2  ♦  £  A0hN 

«  <£  vhNi  - (1  -f  vhN2  +  hNa  +  fvn  • 

But  since  N  »  N(y,&), 

ni  < (1  -£  xi)N  +  r N 
*2  >  <FX1>N  "  R  N 
s3<aoFs+5N  • 

Substituting  into  (5.71)  gives 

|e^(  1/N)  |  <  <£  x^hd  -  £  X,)N  +  (£  X^  i  hN 

-  <’  “Fxi>h<F  V*  +  (1  ~  F  xi  *  R  hN 
+  2  F  V*  +  R  hN  +  F  AohN 

<  2  ^  hN  *  3  £  Anh»  • 

m  n  u 

tj 

Now  since  hN  *  — ,  conclude  that 
n 

l¥1'1"l<tJS{2ff+3jVM))  . 

Similarly, 

|Et(mN+1,<m+1)N)|  <  tj  ;  {2  J  +  3  |  A^mN+1,  (m+1)N)> 

Therefore , 

R  +  3  ^  A^fmN,  (m+1  )N)}tj  ~ 

m“  1 

«  {2  +  3  £  hL(  ^ 

<  t2  +  3  F  A)tJ  R 

<  Utj  , 
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and  let 


V^'i2)  H  *4(tj2)  -  xt(tj1)  -  ;t3* 


By  the  definition  of  approximate  characteristics, 

J2 


where 


-  l  E.(j)  , 

j-jl 


I 

For  j  e  [1,N], 


j  Define 

I 


and  set 
(5.70) 

Then  by  (5.70) 


f(fvh'  if  ‘j" 

I  (|  -  1 )h,  if  aj  > 


|Aj  -  Aj  <  Al(1,N)  B  A0 

»,  ■  (j  <H  !  4j  <|U  , 

N2  -  {j  <  N  s  aj  >  £  Aj}  , 

X-,  “  (0,  1  -  |  A,  -  A0]  , 

X2  “  t1  -  K  X1  +  *0'  ^  ' 

2a  -  n  -fxi  •  V  1  -&  A1  ♦*<>> 

Nn  =  N(Xn,1,N)  ,  n  -  1,2,3 

N,  <  N,  +  N3  , 


Ajjltldt 


Now  we  can  estimate 


sr|w  crl*- 


and  by  definition  of  (J, 


By  (5.51), 


(5.66) 


N  >  N(M,a)  . 


I  -  A,  [t  ,-]  |  <  A  <1,j) 


Thus  for  some  inter ger  m,  0  <  m  <  M-1,  we  must  have 


(5.67) 


A 

Ai(mN,(m+1)N)  <  ^ 


But  this  implies  that  for  some  j  e  [mN,(m+1)N], 


or  equivalently. 


+  Sjh  e  (xA  +  \L(j)k,  xt  +  Al( j )k  +  — )  , 


aj  e  lA^it),  Ajjtt  ]  +  £  , 


where  i  *4(j).  This  follows  because  (5.66)  and  (5.67)  imply  that  there  is  a  fixed  open 
interval  1,  I  c  [A  <j),A  (j)  +  j?')  j-j-  for  all  j  e  [mN,(m+1)N],  such  that  |l|  ' 

where  we  use  (5.62).  Thus  by  Lemma  (5.11), 


and  so 


t)(l,mN,  (nrt-1  )N)  >  {|l|  -  £}  N  >  0  , 


N(I,mN,(m+1)N)  >  1  . 


By  definition  of  approximate  characteristics,  this  implies  that 

iy«Jtyj>  -y»Jvi)| 

for  that  value  of  j  6  [mN,(m+1)N]  for  which  8j  6  !•  Hjus  by  (5.65), 

ivy*1  <3Gos<v  • 

This  verifies  (5.59)  since  tj  <  v  <  tj  implies  | |  ”  |YJl(t;j  +  )|  . 

We  now  verify  (5.60).  Fix  h  <  (5 ( u ) t _  ,  and  set  N  =  tjM  3h  3  so  that 


and  by  (5.64) 


h  “  MN  ' 


N  >  N(M,ft) 


Without  loss  of  generality,  we  do  the  case  l  €  M(tj)\N(tj),  the  case  i  €  W(tj) 
being  similar.  Thus  let  l  e  M(J)\N(tj)  be  fixed.  Define 


•v  v 


then 


l*i<t>|  <  Utj  , 

for  all  i  e  M. 

(Here  (5.59)  says  that  the  strength  of  rarefaction  characteristics  tends  to  zero  uniformly 
with  h  due  to  the  splitting  of  characteristics )  and  (5.60)  says  that  as  h  ♦  0, 
characteristics  move  wih  characteristic  speed. ) 


Proof ; 
(5.61) 


First,  for  l  e  M,  define  [cf  (5.51),  lemma  (5.10)] 


A£(  j  1,  j2) 


j2 

1  G{B(j)+B(j) 

j-jl 


+  C 


°4(j),j} 


Set  A  =  Max{2,GQV0},  so  that 

(5.62)  A4(j1,j2)  <  G0VQ  <  A,  A  >  2  , 

for  all  jl,  j2  <  J.  Let  g  and  y  >  0  be  given.  Let  H  be  the  smallest  integer  such 
that 

(5.63)  M  >  3GQA(2  +  i)y_1  , 


and  define 


<5’64)  “  MN(M,a)  ' 

where  in  addition  to  all  previous  estimates,  Gq  satisfies 

(5.65)  |Y4I  <  G0|Xp(u^[t]  “  Xp(UgCt]  )  |  , 

for  all  l  e  Up,  p  =  1,...,n,  t^  <  t  <  tj  <  tj.  (Recall  that  the  strength  | ( t ) |  is 

ii  0  1 

constant  and  equal  to  |y4|  for  all  t^  <  t  <  t^.)  Such  a  G0  exists  by  Property  (4.3) 
together  with  the  assumption  of  genuine  nonlinearity.  Here  N(M,a)  is  defined  in  Lemma 
(5.11).  We  first  verify  (5.59).  For  this  case  choose  v  >  0,  and  let  t  6  satisfy 
t^  >  v.  Fix  h  <  j$(y)v,  and  set 

N  =  vM-1h-1  , 


so  that 


Theorem  (6.3)  is  a  direct  consequence  of  Proposition  (6.6)  together  with 

PROPOSITION  (6.7);  If  Q(0)  -  Q(T)  <  — - — -  where  u*1  satisfies  our  previous  assumptions 

(GM)2 

then  either  r  e  A/(T)  or  s  e  W(T). 


Proof  of  Theorem  (6.3):  If  Q(0)  -  Q(T)  <  - then  by  Proposition  (6.7)  either 

(GM) 

r  e  W(T)  or  s  e  W(T).  But  by  Proposition  (6.6)  both  r  and  s  are  in  M\V(T).  Thus 

by  contradiction  we  must  have  Q(0)  -  Q(T)  >  — r. 

(GM) 

It  remains  to  give  a  proof  of  Proposition  (6.7).  Proposition  (6.7)  is  a  consequence 

of  the  following  lemmas.  The  idea  is  to  show  that  if  fi(0)  -  ft(T)  <  — ~2»  then  for 

(GM)2 

t  <  T,  Xr  [t]  and  Xs  [t]  are  sufficiently  close  to  X^.  [ 0 )  and  X  [0] ,  respectively,  to 
guarantee  that  the  characteristics  Yr  and  Ya  must  intersect  before  time  T.  Then  by 
Proposition  (5.2),  r  e  N(T)  or  s  e  N(T). 


LEMMA  (6.8) .  Assume  that  Q(0)  -  Q(T)  <  - —.  Then 

(GM) 

(6.70)  l  |y4I  >  |  , 

Ro 

(6.71)  l  |y4!  >£  • 


Proof :  By  Property  (4.2), 

TV{Zp(u£(0))}  -  l  |zp(u*[0])  -  zp(u^[0])|  . 

M 

This  together  with  (6.40)  -  (6.43)  and  (6.46)  implies  that 


\  1**1  • 

rl 

R  L 

(Here  we  use  the  fact  that  z  (u. [t] )  -  z  (u.(t))  is  positive  [negative]  for 

P  *  P  * 

i  e  M+[i  e  M  J,  respectively.)  But 


l  |Ttl 

oL 


l  |7j 

r>L 


1  I 

N(T) 


and  by  Corollary  (5.8) 


since  we  assume  Q(0)  -  Q(T)  <  — — Thus  by  (6.22),  6  <  4-r  ,  so 

(GM) 

rl 

0 

Similarly,  in  the  case  (6.49),  (6.50)  of  Proposition  (6.6),  we  must  have 

l  l*tl  “  l  1**1  "  l  lYtl 

S*  S*  N(T) 

>1.  i  £  >  1 

8M  "  M  M  M  ' 

Since  s  e  Sj,  the  case  (6.47),  (6.48)  immediately  gives  the  conclusion  (6.71).  This 
completes  the  proof  of  Lemma  (6.8). 


LEMMA  (6.9A)  ;  If  1  8  Mp,  S  AKT) ,  p’  j<  p,  and  x^O)  e  (xr(0),  xs(0)),  then  in  time 
[0,T]  Tj  either  intersects  all  the  characteristics  in  Rg  (the  case  p’  <  p)  or  else  it 
intersects  all  the  characteristics  in  S ®  (the  case  p'  >  p) . 

LEMMA  (6.9B);  If  l  eMp,  \  WT),  p'  ft  p  and  intersects  Yr  or  fg  in  time 

T  T 

[0,  — ] ,  then  in  time  [0,T],  Yj  either  intersects  all  the  characteristics  in  RJ  or  else 

it  intersects  all  the  characteristics  in  Sq. 

Proof ;  Since  a  is  assumed  to  be  best  equidistributed,  statement  (5.60)  of  Theorem  (5.12) 
implies  that  for  any  l  eMp\N(T),  1  <  p  <  n,  and  t  <  T, 


(6.72) 


|xt(t)  -  x4(0)  -  Jg  Xa  [t]  dt  I  <  utj  =  (Gghtj/2 


Ix^t)  -  x,(0)  -  X  t|  <  (G^)^  +i-=E  , 


where  without  loss  of  generality  we  have  taken  G0  large  enough  so  that 


-  y 


for  all  l  eMp(J),  t  <  tj. 


v  /  /, 


1 


For  the  proof  of  Lemma  (6.9A),  assume  that  i  eMp A  N(T),  p'  ^  p.  We  do  the  case 
p‘  <  p  and  Yj  intersects  Yr  in  time  [0,  T/2)»  i.e.,  we  show  that  Xj(T)  “  *riiT)  <  0 
for  all  r*  e  [cf.  Prop.  (5.1)].  By  (6.72), 


x  <T)  -  X  (T)  -  X  (T)  -  X  (T)  +  x  (T)  -  x  (T) 


(6.73) 


<  1x^(0)  -  xr(0)]  +  [Xpl  -  Xp]T  +  2E 


+  [xr(0)  -  xr,(0)]  +  [X  -  X  IT  +  2E 


Since  Y4  intersects  Yr  in  time  (0,  T/21,  (6.72)  also  implies 

0  >  x4(|)  -  xr(|)  >  [x4(0)  -  xr(0)J  +  [Xp,  (0)  -  Xp(0)J  |  -  E  . 
Moreover,  by  Lemma  (6.5), 

[xr(0)  -  xr,(0)]  <  |xB  -  xA|  <  y*  . 

Therefore  substituting  into  (6.73)  yields 


x . (T)  -  x  ,  (T)  <  [X  ,(0)  -  X  (0)]  1  +  +  5E} 

Jc  r  p  p  2  o 


,  T  r  6M  V>  5G0T1 

4  -x  2 +  It  +  5<G0hV  2  +  — } 


(6.74) 


{-1  +  **  +  .1°  (G  ht  ^2  +  . 

1  «TX  TX  0  J  MX  '  2 


(6.75) 


By  (6.28)  and  (6.34A)  we  have 

2  EM 


,  2  .  2  ,  1 

6tX  -_2„,  _2  <  3  ' 


6G  MX 


®G2(-t~)X 


by  (6.33),  h  <  (- 


30/2  G‘ 


7l 


)  . 


(6.76) 


10  V  10^  Go2  1/2  1 

ix  <W  2  <  (-557-  )h  <  3 


and  by  (6.34A)  again, 
(6.77) 

Thus  (6.74)  implies 


,0G0  <  _ 
MX  <  3 
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x^tT)  -  xr,(T)  <  0  • 

This  completes  the  proof  of  Lemma  (6.9a). 

For  the  proof  of  Lemma  (6.9B),  assume  l  e  Up,\  N(T) ,  p  p' .  Vie  do  the  case  p'  <  pi 
i.e.,  we  show  that  if  A  e  Upi  *  p'  <  p,  and  x^(0)  e  (xr(0),x),  then  Y^  intersects 
Yr,  in  time  ( 0,T ]  for  ali  r'  e  Rq.  it  suffices  to  show  that  x^fT)  -  xr. (T)  <  0. 

xe(T)  -  xr,(T)  <  [xt(0)  -  xr,(0)J  +  [X^  -  Xp)T  +  2E 
<  ji-  -  XT  +  2E 


-  H  +  Sr +  k  (GohTa>/2  +  srJTX  < 0 


where  we  have  applied  (6. 75)  -  (6.77).  This  completes  the  proof  of  Lemma  (6.9B). 


LEMMA  (6.10). 
(6.69)) 

(6.78) 


Assume  that  Q(0)  -  Q(T)  <  - j.  Then  for  q 

(GM> 


0,  r  and  s  (cf.  (6.63)  - 


Proof .  Write 

By  Corollary  (5.8), 
(6.79) 


since  we  assume  that  Q(0)  -  C(T)  <  - =••  On  the  other  hand,  if  i  e  B_\N(T),  then 

(GM) 

Lemma  (6.9)  implies  that  Y^  intersects  all  the  characteristics  in  either  Rq  or  Sq  in 
time  [Q,T] .  For  example,  assume  the  case  Rq,  and  define 

A*  =  {<A,r'>  ?  l  e  Bq\W(T)  and  r*  e  Rq}  . 

Then  A  #  ±  A  (0)  \  A(T) .  Thus  Lemma  (5.7)  implies 

G. 

(6-80)  l  | Y£)  J Y_,  |  <  ~2~2  * 

A'  gV 
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since  we  assume  Q(0)  -  Q(T)  < 


1 


But  by  Lemma  (6.8) 


(GM) 


2* 


l  IyaI  Ivl  -7  {  i  KIUHy.I} 

i'  *  r  ^  R  v  inn»\  x  _r  x 


(6.81) 


8  \M(T) 


>2M„  X  W*l  • 
8  \V(T) 

<1 


Combining  (6.80)  and  (6.81)  gives 

(6.82) 


1  W.|  < 


2G„ 


*81  2 

8  \N(T)  6G*M 

q 


and  combining  (6.80)  and  (6.82)  gives 

D  *  4G°  <  5 

q  6G2M  M 

where  we  have  applied  (6.22)  and  (6.23).  This  completes  the  proof  of  Lemma  (6.10). 


LEMMA  (6.11)  i  If  Q(0)  -  Q(T)  <  - then 

(GM)* 


(6.83) 

for  q  **  L  and  R. 


Bsq<$  * 


Proof.  First  assume  that  (6.48)  of  Proposition  (6.6)  holds,  and  write 


K  =  l  Iy.I  -  l 


l 


~s  **  » *  i  • 

S'  S"nM(T)  8  ~\N(T) 


By  Corollary  (5.8), 


l  1**1  «<£. 


BenN(T) 

and  moreover  (as  in  the  proof  of  Lemma  (6.10)),  A'  £.  A(0)\A(T)  where 

A'  =  {<f,s'>  !  *  e  B”  ,  s'  e  Sg,  xs, (0 )  *  xs<0)}  . 

(Here  we  apply  Lemma  (5.2).)  Thus  as  in  (6.80)  -  (6.82),  Lemma  (5.7)  implies  that 

4G 

‘  8- '  ^  ‘  ’  ""B  • 


-56- 


where  we  apply  (6.48)  in  place  of  Lemma  (6.8). 

Now  assume  (6.49)  -  (6.52)  of  Proposition  (6.6)  holds.  In  this  case  (6.50)  -  (6.52) 

f 

satisfy  the  hypotheses  (5.42),  (5.43)  of  Proposition  (5.9)  with  L  «  jj.  Since  we  assume 


Q(0)  -  Q(T)  <  — 1— , 


1  1  2  6  2  6  2 
- j,  and  - —  *  6  (— )  <  G  (— )  >  we  can  conclude 

(GM)  (GM)2  "  0  M. 

8-  <  A 

sq  M 


for  q  »  L, R.  This  completes  the  proof  of  Lemma  (6.11). 


LEMMA  (6.12)8  Assume  that  Q(0)  -  Q(T)  <  - 5.  Then 

(GM)2 

(6.84)  X  [0]  -  X  (0]  >  -i-  . 

r  s  GQM 


Proof.  We  prove  Lemma  (6.12)  by  satisfying  the  hypotheses  of  Lemma  (6.2).  First,  (6.46) 


implies 

(6.85) 

and  (6.47),  (6.49)  imply 


(6.86) 


Zp(u^Ol)  , 


Zp(u£[0})  <  Ijj  . 


Let  u  =  u£[0] ,  and  let  S  denote  the  p-shock 


'i  0  * 

s 


Let  {u  ,uR}  denote  the  left  and  right  states,  and  a  the  speed  of  the  p-shock  S.  By 

(5.49)  X^ [0]  •  a,  and  (6.86)  gives 

0  <  z  <uR>  <  z  (uR[0] )  <  i-  . 
p  p  s  8M 


By  (6.39)  we  must  have 


VuL)  4  i  • 


z  (u11)  +  z  (uR) 


By  (6.85) 


z  (uL)  +  z  (uR)  , 

a  (u)  -  -£ _ P  -  >  —L 

pl  2  16M  • 
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Moreover,  for  p‘  f  p,  the  difference  in  zp,  between  uR  end  u  ie  bounded  by  the 
total  variation  in  zp,  of  all  wavee  that  lie  between  xr(0)  and  x#(0)  at  t  •  0.  By 
Property  (4.3),  this  can  be  estimated  by 

l*p.<«>  -  *p.<uR)|  *  B0  +  1  |zp,(u*[01)  -  *p,  <Uj[0)  )  | 

M 

P 

M  M  M  ' 

where  we  have  applied  Lemma  (6.4)  and  Lemma  (6.10)  in  the  second  inequality.  Thus  Lemma 
(6.2)  applies  with  L  *  -j  to  giws 

*  <u>  ~  o  >  ~  • 
p  GqM 

Therefore  we  conclude 

*10]  -  A  [0]  -  A  (u>  -  o  >  ~  . 

r  s  p  G0M 

Proof  of  Proposition  (6.7):  Assume  that  Q(T)  -  Q(0)  <  — — =-.  He  show  that  Y 

(GM) 

intersects  Tg  in  time  (0,  T/2] .  By  Proposition  (5.2)  this  implies  that  either  r  e  N(T) 
or  s  e  N(T).  Thus  it  suffices  to  show  that  xb(t/2)  -  XjlT/2)  <  0.  By  (5.53), 

|*alT/2]  -  Ag[0]  |  <  Gq{Bb  +  B"  +  [Q(0)-Q(T)]1/2  } 


(6.87) 


(6  6  1  i  3G0® 

*  G0^M  +  M  +  GM^  *  M  ' 


where  we  applied  Lemmas  (6.10)  and  (6.11).  Therefore  we  can  use  (5.42)  and  (5.44)  to 


obtain 


(6.88) 


(T/2)  <  xa<0)  +  f£/2  Ag(t)dt  +  (G0ht7)1/2 


T  0  T  Vo 

<  V0)  +  xsl°l  7  +  TT2-  J+  (Gohtj>  2  • 


Similarly, 


(6.89) 


T  3®G0  T  Vo 

Xj. (T/2 )  >  *r(0)  +  Ar[0]  ^  ?  -  (G0htj/2  . 


Subtracting  (6.89)  from  (6.88)  gives 


x  (T/2)  -  x  (T/2)  <  [x  (0)  -  x  (0)]  +  {A  [0]  -  A  [0]}  ± 
s  r  s  r  s  r  2 


r  3OG0T  Vo  , 

+  "Hm-  +  (G0hV *  1 


>>>>•/ 


•«  ->  v\.\h  *  -  *  .  V  •  .  Vs  .  ».*  .* 


*  *  -*  A*.**".*,f.  < 

*  .  •  •  .  *  v"  *,***•••*.  * 

■i-.' 


v  V 


and  by  Leninas  (6.5)  and  (6.12)  respectively. 


|x  (0>  -  x  (0||  <  |S  . 


V",  -  V01  '  -  sjs 


Therefore 


3fiG  T  1 

i(T/2)  -  Xr(T/2)  <  {“  +  -±-  +  2(G0htJ)/2  - 


(6.91) 


[V*  ,  [r2.  .  Wa’1 

*  5t  +  “o4  + - T - 


By  (6.28)  and  (6.22),  (6.23), 


(6.92) 


By  (6.22) 


(6.93) 


2G0M2e  2G0M2£  2Gq  3  , 

“^“'^V^“26Go<2  * 


6G26  <  1  . 


By  (6.28)  and  (6.33), 


(6.94) 


4G.M(G  ht  //2  4/2  m3/2fr  4/2  MG3/2  V, 

0  0  J  .  0  f  0  w  T* 


<  < - o  i  .  i 

•Z2  12/2  mg3/2  3 


Finally,  putting  (6.92)  -  (6.94)  into  (6.91)  yields 

xg(T/2)  -  Xj, (T/2)  <  0  . 

This  completes  the  proof  of  Proposition  (6.7),  and  hence  alBO  completes  the  proof  of 

Theorem  (6.3).  The  proof  of  Theorem  (6.3)  also  applies  to  the  case  of  periodic  data,  in 

which  case  I  I  .,  I  I  and  Q(t)  are  defined  on  each  period. 

L  S 
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